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1.0 SUMMARY 




I 

I 

I 


This User's Manual contains a complete description of the computer codes 
developed and utilized under NASA contract NAS3-23716 "Assessment of a Three 
Dimensional Boundary Layer Code to Predict Heat Transfer and Flow Losses in a 
Turbine". The SETUP code calculates a general nonorthogonal surface coordinate 
system and the boundary layer edge conditions from a known static surface pres- 
sure distribution. The BL3D code utilizes this coordinate system and the edge 
conditions for the calculation of the three dimensional boundary layer on an 
arbitrary surface such as a turbine blade or end wall. A companion report to 
this User's Manual which is entitled "Assessment of a 3-D Boundary Layer Analysis 
to Predict Heat Transfer and Flow Field in a Turbine" documents the technical 
approach used in the SETUP and BL3D codes and presents results which have been 
obtained in the application of the three dimensional boundary layer theory to the 
flow in a turbine passage with blade rotation. Although the primary thrust of 
this work has been for gas turbine application, the compressible three dimen- 
sional boundary layer analysis is very general and can be used for a wide range 
of applications in internal and external aerodynamics. 

The User's Manual is essentially divided into two parts. In the first part, 
Sections 3.0, A.O, and 5.0, a description of the input, output, and general oper- 
ation of both codes is presented. In the second part. Sections 6.0 and 7.0, a 
detailed description of the computer codes including a listing and description of 
all FORTRAN variables in COMMON blocks and a detailed description of the subrou- 
tines are given. These latter sections of the User's Manual are presented for 
the user who wishes to obtain a detailed understanding of the codes. 
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2 . 0 INTRODUCTION 

This User's Manual describes the computer codes known collectively as the 
Three Dimensional Boundary Layer (BL3D) code wich was developed under NASA 
Contract NAS3-23716. A description of the analysis used in the SETUP code is 
given by Anderson (Ref. 1). A comprehensive outline of the analysis used in the 
BL3D code is also given by Anderson (Ref. 1). The three dimensional boundary 
layer analysis used in the BL3D codes was developed by Vatsa (Refs. 2 and 3). 
Reference 3, in particular, describes in detail the solution algorithm including 
the finite difference equations and the matrix inversion method. Assessment of 
these codes for gas turbine applications is also given in Ref. 1 through 3. The 
SETUP code contains a prepackaged surface spline fitting routine developed by 
McCartin (Ref. 6) which is not described. 

This report is divided into sections for the convenience of the user. 

Section 3.0 contains an overall description of the computer codes describing how 
they are linked together functionally and through data files. Sections 4.0 and 
5.0 contain a description of the input and output files so that the user can run 
the codes and interpret the printed output. These sections also contain a sample 
input data file to assist the user in setting up new cases. Although these codes 
are very robust, they sometimes fail. If the cause of the failure is known, a 
DIAGNOSTIC is printed and the calculation is terminated. A description of these 
failures and DIAGNOSTICS is also given. Sections 3.0 through 5.0 contain all the 
information required by the user who wishes to run these codes but does not 
require a detailed understanding of the codes. 
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Sections 6.0 and 7.0 contain a more detailed description of these computer 
codes for the user who wishes to modify these codes. These sections are divided 
into three subsections. The first subsection describes the MAIN program 
including a global tree structure chart which provides the user with an overview 
of the code. This section includes all calling and called subroutines. The 
second subsection contains a list of all FORTRAN symbols which are contained in 
COMMON blocks. These FORTRAN symbols may be identified with the corresponding 
algebraic symbols which are the same as those used in Ref. 1. The third 
subsection contains a detailed description of each subroutine. These subroutine 
descriptions are arranged according to the same format. First, the object of the 
subroutine is described. Second, any input options used in the subroutine are 
described so that the user can understand any branching that takes place. Third, 
a list of local FORTRAN symbols is given. Finally a description of the analysis 
used in the subroutine is given. If the analysis is described in Ref. 1, it is 
referenced by subsection. If the analysis is not described in Ref. 1, it is 
described in these subroutine descriptions. In cases where well known analyses, 
such as interpolation formulas, are used the source reference is given. 
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3.0 GENERAL DESCRIPTION OF BL3D CODE 

The Three dimensional boundary layer (BL3D) code requires three separate 
analyses: 1) an analysis to construct a general nonorthogonal coordinate system 

for a twisted turbine blade, 2) an analysis referred to as a surface Euler 
analysis to calculate the boundary layer edge conditions from a known wall static 
pressure distribution, and 3) an analysis to calculate the three dimensional 
boundary layer growth on the turbine blade surface. These analyses are contained 
in two separate computer codes. The SETUP code contains the coordinate analysis 
and the boundary layer code analysis. The BL3D code contains the boundary layer 
analysis. A detailed description of the coordinate analysis and boundary layer 
edge analysis contained in the SETUP code is given by Anderson (Ref. 1). A 
general outline of the boundary layer analysis contained in the BL3D code is 
given by Anderson (Ref. 1). A detailed description of the boundary layer 
analysis including the finite difference equations and the matrix inversion 
analysis is given by Vatsa (Refs. 2 and 3). 

A flow chart describing the operation of these computer codes is shown on 
Fig. 1. The geometry analysis calculates a surface coordinate system which is 
required by both the surface Euler analysis and the three dimensional boundary 
layer analysis. The surface Euler analysis calculates the boundary layer edge 
conditions. This data is then input into the 3D boundary layer analysis. A 
brief description of these codes is given below. 
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Geometry Analysis 

Since the boundary layers lie on the blade surface, a convenient coordinate 
system is a surface coordinate system (X1,X2,X3) which has the properties that 
two coordinates (X1,X2) lie on the blade surface and the third coordinate (X3) is 
normal to the blade surface. XI is generally taken to be in the streamwise 
direction and X2 is generally in the crossflow direction. This coordinate 
system, in general, must be a non-orthogonal coordinate system, therefore XI is 
not necessarily orthogonal to X2 . Furthermore we note that a given twisted 
turbine blade is not a developable surface (i.e., it can not be flattened to a 
plane surface without distortion). Thus if we describe the turbine blade in 
Cartesian coordinates (Y1,Y2,Y3), where Y1 is generally in the streamwise 
direction, Y2 in the crossflow direction, and Y3 orthogonal to (Y1,Y3), then a 
mathematical transformation must be found for the relation between the computa- 
tional coordinates (X1,X2,X3) and the physical coordinates (Y1,Y2,Y3). In 
general, the equation of the physical surface can be expressed as Y3 = Y3(Y1,Y2). 
Thus if a unique relation can be found Y1 = Y1(X1,X2) and Y2 = Y2(X1,X2), then 
the third coordinate is known Y3 = Y3(X1,X2). This unique relation is provided 
by the transfinite mapping of Gordon and Thiel (Ref. 4). Once the transformation 
from the computational coordinates (X1,X2,X3) to the physical (Cartesian) 
coordinates is known, the covariant matrix tensor, metric scale coefficients, and 
the direction cosines from the computational coordinates to the Cartesian 
coordinates can be determined using the relations given in Ref. 5. In addition 
it is noted that since the turbine blade surface is only known as numerical data, 
the surface spline analysis of McCartin (Ref. 6) and Spath (Ref. 7) is used to 
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interpolate points on the surface and calculate first and second derivatives of 
the surface. Details of this geometric analysis are given by Anderson 
(Ref. 1). 

Surface Euler Analysis 

The three dimensional boundary layer equations require the boundary layer 
edge conditions, which include the two free stream velocity components, their 
derivatives, and the free stream temperature, as boundary conditions over the 
whole computational domain. These may be obtained directly from experimental 
data, however, in practice it is very expensive to survey an entire three 
dimensional flow field. These conditions may also be obtained from a solution of 
the Euler equations. However, since the generation of the passage vortex is 
viscous in nature, the Euler solutions may not provide good edge conditions for 
the crossflow without special treatment. A third method is to obtain the edge 
conditions from the experimental pressure distribution by intergrating the 
surface Euler equations. These surface Euler equations are obtained from the 
boundary layer equations in the limit as the normal coordinate approaches the 
edge of the boundary layer so that normal derivatives of all variables go to 
zero. Since the pressure distribution is known, only two momentum equations and 
an energy equation are required for a complete solution of all the dependent 
variables. These surface Euler equations are hyperbolic in character and may be 
solved using the same computational (X1,X2) grid, finite difference module, and 
inflow edge conditions as the three-dimensional boundary layer analysis. A 
detailed presentation of this analysis are given by Anderson (Ref. 1). 
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3-D Boundary Layer Analysis 

The three dimensional boundary layer analysis used in this computer code is 
that developed by Vatsa (Refs. 2 and 3). This analysis solves the finite 
difference form of the compressible three-dimensional boundary layer equations 
(including the energy equation) in a non-orthogonal surface coordinate system 
which includes Coriolis forces produced by coordinate rotation. These equations 
are solved in Levy-Lees variables using an efficient, implicit, fully coupled 
finite differnce procedure. These boundary layer equations are hyperbolic in 
character and require, besides the edge boundary conditions, the inflow condi- 
tions along the computational boundaries where flow enters the domain. These 
inflow conditions are estimated by solving the 3-DBL equations along the inflow 
boundary using an assumption of either local similarity or plane of symmetry 
which reduces the equations to two dimensional differential equations with edge 
conditions provided by the boundary layer edge analysis. 

Operation of Computer Codes 

The description of the input and output of the SETUP code is given in 
Section 4.0 and a description of the input and output for the BL3D code is given 
in Section 5.0. Sample inputs for both codes are given in the respective 
sections. Since the BL3D code requires the surface coordinate data and the 
boundary layer edge conditions, the SETUP code must be run first. The linking of 
these two codes through data files is shown in Fig. 2. Note on Fig. 2 that 
output Units 8, 9, and 10 for the SETUP code correspond to the input Units 10, 

11, and 12 respectively for the BL3D code. Output Units 13, 15, and 16 are 
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special output files for CALCOMP or TEKTRONIX plot packages. Since plotting 
codes are machine dependent, they are not available on the NASA Lewis version of 
these codes. Unit 20 for the BL3D code serves as both an input and output file. 
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Fig. 1 Flow Chart of Geometry, Surface Euler and 3-D Boundary Layer Analysis 
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4.0 OPERATION OF SETUP CODE 
4.1 Input Format 

The input format is described on the pages which follow. The type and 
amount of data read by the code is controlled by the input options lOPTl through 
I0PT12 and the user should be very careful to be sure that the input data agrees 
with the input options. In general the input cards are loaded as follows: 


Card( s) 

0 

Title Card 


Card( s) 

1 

Computational Grid Data 


Card( s) 

2 

Computational Boundary Data 


Card( s) 

3 

Coordinate Surface Data 


Card( s) 

4 

Pressure Distribution Data 


Card( s) 

5 

Inflow Conditions Data 


Card 

6 

Print Options 


Card( s) 

7 

Reference Flow Variables 


Card( s) 

8 

Rotor Parameters (load only 

if I0PT8 = 1) 

Card( s) 

9 

Upstream Flow Conditions (load 

only if lOPTll = 1) 


In the loading of the input cards, the variable name, column numbers, and 
description are given. Integer numbers are right adjusted and real numbers may 
have any number of significant digits. 
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Card(s) 0 Title Card (load 72 alphanumeric characters) 

Card(s) 1 Computational Grid Data 

^a_rd_ _ 

I0PT4 1-10 Mesh distribution option 

I0PT4 = 0 Calc, uniform grid 

1 Calc, stretched grid 

2 Read grid points 

LPLAST 11-20 No. XI grid points 

KPLAST 21-30 No. X2 grid points 

^a^d_ _ Ji.b_Sjt r^t£^h^n^ ^a^ame^er^ 

A 1-10 Stretching parameter 

^ar^d^s^^ J^b_Ijip_ut_g_ri^ P^oan^s 

(Xl(LP) ,LP=1 , LPLAST) XI grid points 

Load LPLAST numbers (8 numbers/card) 

(X2(KP) ,KP=1 , KPLAST) X2 grid points 

Load KPLAST numbers (8 numbers/card) 


Card(s) 2 Computational Boundary Data 


^a^d_ _ ^a_0^tj^o£s 
I0PT2 1-10 


LPLAST 11-20 
KPLAST 21-30 


(8110) 


Grid boundary option 
I0PT2 = 0 Cartesian Grid on a plane 
= 1 Skewed grid on a plane 
= 2 ID mapping Y2=Y2(X1,X2) 

= 3 ID mapping Y1=Y1(X1,X2) 

= 4 2D mapping 

No. XI input boundary points 
No. X2 input boundary points 


IV- 2 


\ 

I 

1 



R85-956834-1 


£a^d_ _ ^b_B^und^r^ 


YIL 

1-10 

Lower 

Y1 

boundary 

YIU 

11-20 

Upper 

Y1 

boundary 

Y2L 

21-30 

Lower 

Y2 

boundary 

Y2U 

31-40 

Upper 

Y2 

boundary 

SSIO 

41-50 

Initial . 

arc length on XI 


£a^d 


YlOO 

1-10 

Y1(0,0) 



YlOl 

11-20 

Y1(0,1) 



Ylll 

21-30 

Yl(l,l) 



YllO 

31-40 

Y1(1,0) 



Card 2c 

Boundary 

Corners (4F10.0) 


_IOPT2_=_ 

Y200 

1-10 

Y2(0,0) 



Y201 

11-20 

Y2(0,l) 



Y211 

21-30 

Y2(l,l) 



Y210 

31-40 

Y2(l,0) 



Card 2b 

Y1 Boundary Points (2F10.0) 


_I0PT2_=_: 

YIL 

1-10 

Lower Y1 boundary 



YIU 

11-20 

Upper Y1 boundary 
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Cavd^l 2c_ ^ide_2_and_Sj.d£ 4 _ 18F10.0) 

(YlXll(LB) ,LB=1,LBLAST) Y1 boundary points side 2 (X2 = 1) 

(Y2X11(LB),LB=1,LBLAST) Y2 boundary points side 2 (X2 = 1) 

(Y1X10(LB),LB=1,LBLAST) Y1 boundary points side 4 (X2 = 0) 

(Y2X10(LB),LB=1,LBLAST) Y2 boundary points side 4 (X2 = 0) 

Load 4*LBLAST numbers (8 numbers /card) 

^ar^d_ _ 2b_ I.2_Bo^undar2 Poin^s_ 12F1^. 0) f. 1) 

YIL 1-10 Lower Y1 boundary 

YIU 11-20 Upper Y1 boundry 

jCa£dU^ ^c_ lile_l_and_Sj.d£ 2 _ ifll°P12_=_3l 

(Y10X2(KB),KB=1,KBLAST) Y1 boundary points side 1 (XI = 0) 

(Y20X2(KB),KB=1,KBLAST) Y2 boundary points side 1 (XI = 0) 

(Y11X2(KB) ,KB=1,KBLAST) Y1 boundary points side 3 (XI = 1) 

(Y21X2(KB),KB=1,KBLAST) Y2 boundary points side 3 (XI = 1) 

Load 4*KBLAST numbers (8 numbers/card) 

2a£d2s2. 2,^_P^ur^ 

(Y10X2(KB),KB=1,KBLAST) Y1 boundary points side 1 (Xl = 0) 

(Y20X2(KB),KB=1,KBLAST) Y2 boundary points side 1 (XI = 0) 

(YlXll(LB) ,LB=1 .LBLAST) Y1 boundary points side 2 (X2 = 1) 

(Y2X11(LB),LB=1,LBLAST) Y2 boundary points side 2 (X2 = 1) 

(Y11X2(KB),KB=1,KBLAST) Y1 boundary points side 3 (XI = 1) 

(Y21X2(KB),KB=1,KBLAST) Y2 boundary points side 3 (Xl = 1) 

(Y1X1((.(LB) ,LB=1 .LBLAST) Y1 boundary points side 4 (X2 = 0) 

(Y2Xl(fi(LB) ,LB=1 jLBLAST) Y2 boundary points side 4 (X2 = 0) 

Load 4*(LBLAST + KPLAST) numbers (8 numbers/card) 
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Card(s) 3 Coordinate Surface Data 

^a^d_ _ ^a_0^t^ojis 

lOPTl 1-10 Surface geometry option 

lOPTl = 0 plain surface 
= 1 Y3 = Y3(Y1) 

= 2 Y3 = Y3(Y2) 

= 3 Y3 = Y3(Y1,Y2) 

LSLAST 11-20 No. input Y1 points 

KSLAST 21-30 No. input Y2 points 

^a^d^s_)^ ^b_S^r^ac^e_Coo^dJ^nat^s_ Z. J. 

(Y1(LS),LS=1, LSLAST) Y1 coordinate 

(Y3(LS),LS=1, LSLAST) Y3 coordinate 

Load 2*LSLAST numbers (8 numbers/card) 

^a^d^s_). ^b_Sjir^a«_Coord^n^te^s_ *.2.^ Z. 1. 

(Y2(KS),KS=1, KSLAST) Y2 coordinate 

(Y3(KS),KS=1, KSLAST) Y3 coordinate 

Load 2*KSLAST numbers (8 numbers/card) 

Cajrd^s_)^ ^b_S^r^ac^e_Co^o^dJ^nate^s_ ^.2 Z. 2. 

(Y1(LS),LS=1, LSLAST) Y1 coordinate 

Load LSLAST numbers (8 numbers/card) 
(Y2(KS),KS=1, KSLAST) Y2 coordinate 

Load KSLAST numbers (8 numbers/card) 
((Y3(LS,KS),LS=1, LSLAST), KS=1, KSLAST) Y3 coordinate 
Load LSLAST*KSLAST numbers (8 numbers/card) 
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CARD(S) 4 Pressure Distribution Data 


£ar^d_ _ ^a_ ^ 50 ^^ 

I0PT3 1-10 Pressure data option 

I0PT3 = 0 No pressure data 
= 1 CP=CP(X1,X2) 

= 2 CP=CP(Y1,Y2) 

= 3 CP=CP(Y1,X2) 

= 4 CP=CP(X1,Y2) 


LILAST 

11-20 

No. XI or Y1 input points 



KILAST 

21-30 

No. X2 or Y2 input points 



I0PT7 

31-40 

SMOOTH CP DATA I0PT7 times 



I0PT8 

41-50 

Rotor/Stator option 
I0PT8 = 0 Stator 

= 1 Rotor load data 
=-l Rotor load data 

hub 

tip 

to tip 
to hub 

Card(s) 

4b Coordinate 

Data (8F10.0) 



(YII(LD) 

,LD=1 .LILAST) 

Y1 Coordinate 



Load LILAST numbers (8 numbers/card) 


(Y2I(KD) ,KD=1 ,KILAST) Y2 Coordinate 

Load KILAST numbers (8 numbers/card) 
((Y3I(LD,KD) ,LD=1 , LILAST) ,KD=1 jKDLAST) CP (Y1,Y2) 

Load LILAST*KILAST numbers (8 numbers/card) 
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Card(s) 5 Inflow Conditions Data 

^a£^d_ _ ^ 

I0PT5 1-10 Inflow option 

I0PT5 = 0 UlE = SQRT(l-CP) : U2E =0.0 
I0PT5 > 0 integrate bl edge equations 

= 1 U2(0,X2) = 0,0 : U2(X1,0) = 0.0 

= 2 U2(0,X2) = 0.0 : U2(X1,0) Input 
= 3 U2(0,X2) Input : U2(X1,0) = 0.0 
= 4 U2(0,X2) Input : U2(X1,0) Input 

LULAST 11-20 No. XI input points 

KULAST 21-30 No. X2 input points 

(YI1U(LP),LP=1, LULAST) Y1 coordinate 

Load LULAST numbers (8 numbers/card) 

(U2X1(LP),LP=1, LULAST) U2(Y1) 

Load LULAST numbers (8 numbers/card) 

Ca£dXs2 lb_X2=0 £nf low_b£undar£ iSm.O) I_f ^OPT^ £ 2 

(YI2U(KP),KP=1, KULAST) Y2 coordinate 

Load KULAST numbers (8 numbers/card) 

(U2X2(KP) ,KP=1, KULAST) U2(Y2) 

Load KULAST numbers (8 numbers/card) 


IV- 7 



R85-956834-1 

^a_r<Ks2^ _5^b_X2^=^ _inf l^v^b^und^r^ ^1. Z. iL 

(YIIU(LP) ,LP=1,LULAST) Y1 coordinate 

Load LULAST numbers (8 numbers/card) 

(U2XKLP) ,LP=1,LULAST) U2(Y1) 

Load LULAST numbers (8 numbers/card) 

(YI2U(KP),KP=1,KULAST) Y2 coordinate 

Load KULAST numbers (8 numbers/card) 

(U2X2(KP) ,KP=1,KULAST) U2(Y2) 8 numbers/card) 

Load KULAST numbers (8 numbers/card) 

Card 6 Print Options (5110) 

I0PT6 1-10 Print every lOPTbth station 

I0PT9 11-20 Coordinate print option 

I0PT9 = 0 Output Y1 and Y2 
= 1 Output SI and S2 

lOPTlO 21-30 Set lOPTlO = 1 

lOPTll 31-40 Upstream conditions 

lOPTll = 0 Uniform conditions 

= 1 Nonuniform conditions 

I0PT12 41-50 Print data on K = I0PT12th line 


IV- 8 


R85-956834-1 


Card(s) 7 

Reference 

Flow Variables (6F10.0) 


PINFI 

1-10 

Static pressure 

(psf) 

TINFI 

11-20 

Static temperature 

(deg. R.) 

UINFI 

21-30 

Velocity 

(ft/sec) 

BETH 

31-40 

Not used 


XLREF 

41-50 

Length 

(ft) 

DHINF 

51-60 

Not used 


Card(s) 8 

Rotor Parameters (Load only if I0PT8 

= 1) 

Card_ 

8a Rotation Vector (4F10.0) 


VBTI 

1-10 

Rotor tip speed 

( ft/sec) 

GOl 

11-20 

Components of omega 


G02 

21-30 

in the Yl ,Y2, Y3 


G03 

31-40 

directions 


Ca£d_ 

8b Rotor 

Parameters (6F10.0) 


RADHI 

1-10 

Hub radius 


RADTI 

11-20 

Tip radius 


YRADHI 

21-30 

Hub Y2 coordinate 


YRADTI 

31-40 

Tip Y2 coordinate 


GAP I 

41-50 

Mid span gap 


DELCP 

51-60 

Not used 
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Card(s) 9 Upstream Flow Conditions (load only if lOPTll - 1) 


Y2RI(K) 

1-10 

Spanwise coordinate 

(Y2) 


PTAI(K) 

11-20 

Inlet 

absolute total 

pressure 

(psf) 

TTAI(K) 

21-30 

Inlet 

absolute total 

temperature 

(deg. R.) 

UZAI(K) 

31-40 

Inlet 

axial velocity 


(ft/sec) 

BTRI(K) 

41-50 

Inlet 

flow angle 


(deg. ) 


Load KILAST cards 
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4.2 Output Foraat 

The output for the SETUP code echoes the input data and also lists certain 
additional calculated data which can be used for input to the BL3D code. The 
output description is given by headings on the following pages. 


Unit 12 Output (file name OUTPUT.DAT) 


Title Card 


72 alphanumeric characters 
Input Options 

All input options are listed in order with its value and meaning beside 
it . 

Independent Variables 

The input calculation grid points X1(L) and X2(K) are listed as input. If 
I0PT4 =0 or 1, the calculated grid points are listed. 

Computational Boundaries 

The input computational boundaries are listed by sides where: 

Side 1 Y1(0.0,X2) , Y2(0.0,X2) 

Side 2 Y1(X1,1.0) , Y2(X1,1.0) 

Side 3 Y1(1.0,X2) , Y2(1.0,X2) 

Side 4 Y1(X1,0.0) , Y2(X1,0.0) 

are parametric pairs of points describing the boundaries. If IOPT2 < 4 , the 
boundaries are calculated and then listed. 
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Coordinate Surface Data 

The Cartesian coordinates of the computational surface are listed next. It 
is assumed that the independent variable Y1 is only a function of L and the 
independent variable Y2 is only a function of K. These arrays are listed first. 
Next Y3(Y1,Y2) is listed. For lOPTl = 0 , the computational surface is a plane 
surface and Y3 is not listed. For lOPTl =1 or 2, the surface is a two dimen- 
sional (developable) surface Y3=Y3(Y1) or Y3=Y3(Y2). With these options only Y1 
and Y3 are listed or Y2 and Y3. 

Pressure Distribution Data 

The independent variables are listed first. Again it is assumed that Y1 or 
XI is only a function of L and Y2 or X2 is only a function of K depending upon 
the input option I0PT3. The pressure coefficient data is listed next. 

Inflow Conditions (If I0PT5 > 1) 

The independent variables Y1 or XI are listed first and followed by the 
inflow velocities U2(Y1) or U2(Xl). Then the independent variable Y2 or X2 is 
listed followed by U2(Y2) or U2(X2) . 

Flow Parameters 

The following input parameters are listed: 


PINFI 

P» 

Reference 

static pressure 

(psf) 

TINFI 

too 

Reference 

static temperature 

(deg. R.) 

UINFI 

Uoo 

Reference 

velocity 

( ft/ sec) 

BETH 

^00 

Reference 

inlet flow angle 

atan(UZA/UPR) 
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XLREF 

£ 

Reference length 

DHINF 


Not used 

VBTI 

Vb 

00 

Rotor tip speed 

GOl ') 

1 

Direction cosines of rotation 

G02 ' 


vector with respect to the 

G03 

) 

1 

Cartesian coordinates 

RADHI 

% 

Hub Radius 

RADTI 

Rj, 

Tip Radius 

YRADHI 

Yr 

Hub Y2 coordinate 

YRADTI 

Yt 

Tip Y2 coordinate 

GAPI 


Reference gap 

DELCP 


Not used 


(ft) 


( ft/ sec) 


The following calculated parameters are listed next. These quantities are 
nondimensionalized in the same manner as the input to the BL3D code and may be 
used as input to the BL3D code. Note that rotor parameters are only listed if 


I0PT8 > 0. 


PINF 

Pco 

Reference 

freestream static pressure 

TREF 

Tref 

Reference 

static temperature 

TINF 


Reference 

freestream static temperature 

HINF 


Reference 

freestream static enthalpy 

DHINF 


Not used 


XMINFI 

00 

Reference 

freestream Mach number 

OMEGA 

Q. 

Magnitude 

of rotation vector 
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OMEGAl \ 

OMEGA2 \ Cartesian components of rotation vector 

0MEGA3 ^ 

VBINF Vg^ Reference rotor speed 

UPINF Reference inlet tangential velocity 

ALPII OjQQ 90. - BETH 

HTABSF H.p„ Reference absolute total enthalpy 

PTABSF Reference absolute total pressure 

Nonuniform Inflow Conditions (if I0PT8 = 1) 

The following inflow conditions are input and listed. 

Y2RI(K) Yj Y2 coordinate 

PTAI(K) Inlet absolute total pressure (psf) 

TTAI(K) t.j. Inlet absolute total temperature (deg. R) 

UZAI(K) Inlet absolute axial velocity (ft/sec) 

BTRI(K) Inlet relative flow angle atan(UZA/UPR) (deg.) 

Rotor Inlet Conditions (if I0PT8 = 1) 

The following variables are calculated and listed. 

Y2RI(K) Y 2 Y2 coordinate 

PT/PINF P'p/Pjoo Absolute total pressure ratio 

TT/TINF T^/T^„ Absolute total temperature ratio 

P/PINF P/P<xj Static pressure ratio 

T/TINF T/To 5 Static temperature ratio 
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VBB 

Vb 

Rotor velocity 


BTA 

^1 

Absolute 

inlet 

flow angle atan(UZA/UPA) (deg.) 

UlA 

/V/ 

Ul 

Absolute 

inlet 

velocity 

UPA 

% 

Absolute 

inlet 

tangential velocity 

UZA 


Absolute 

inlet 

axial velocity 

UlR 

Ul 

Relative 

inlet 

velocity 

UPR 

% 

Relative 

inlet 

tangental velocity 

QMS 

«s 

Relative 

inlet 

vorticity 

U2S 

^2s 

Estimated 

inlet 

; spanwise velocity 


Boundary Layer Edge Conditions 

The boundary layer coordinates and edge conditions are calculated on a 
LPLAST*KPLAST grid. Depending on the input option I0PT6, the following quantities 
are listed at grid points which are multiples of I0PT6 . 


L , K 


Coordinate point 

XIP , X2P 

Xi. X^ 

Computational coordinates 

YIP , Y2P 

Yi. Y^ 

Cartesian coordinates 

CP 

S 

Pressure coefficients 

UlE , U2E 

Ule> ^2e 

Boundary layer edge velocity components 

DUIEDXI 

au^g/ 8Xj ^ 


DU1EDX2 

auig/ax2 


DU2EDX1 

SDje/aXi 1 

^ Boundary layer edge derivatives of the 
i velocity components 

DU2EDX2 


1 
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0MEGA3 

% 

Vorticity component nomal to surface 

P/PINF 

P 

Static pressure ratio 

H/HINF 

I 

Total enthalpy (rothalpy) ratio 


Unit 11 Output (file name CORSCL.DAT) 

The output on unit 11 provides a detailed short listing of all L stations 
for a given K = KPT station. This output is very useful for aligining a computa- 
tional coordinate with the corresponding Cartesian coordinate and the arc length 
along the coordinate surface. It also lists the coordinates of a calculated 
streamline starting at the point (L = 1 , K = KPT). The pair of points (S1P,SLS) 
form a parametric curve of the streamline on the surface starting at the point (L 


= l.K = 

KPT) . The 

following variables are currently listed on the file. 

L 


Coordinate point 

XIP 

Xl 

Computational coordinate 

YIP 

^1 

Cartesian coordinate 

SIP 

Si 

Arc length along XI coordinate 

Y3P 

^3 

Cartesian coordinate of surface 

CP 

Cp 

Static pressure coefficient 

UlE 

Ule 

Edge velocity components 

U2E 

«2e 


0MEGA3 

% 

Vorticity component normal to surface 

SLS 

S2s 

Spanwise arc length of streamline 

ALPE 


Flow angle at streamline 
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4.3 Output Data Files 

The output data files are written on units 8, 9, and 10. These data files 
contain some of the data which is required as input by the BL3D code. A detailed 
listing of the files are given below. 

Unit 8 Metric Data (7E14.7) 

HHl Metric scale coefficient for XI 

HH2 H 2 Metric scale coefficient for X2 

DHll 9Hj/3Xj^ D(HH1)/DX1 

DH12 3H^/3X2 D(HH1)/DX2 

DH21 3H2/3Xj D(HH2)/DX1 

DH22 3H2/3X2 D(HH2)/DX2 

GL12 Gj ^2 Covariant metric tensor component 

DG |G| Determinant of covariant metric tensor 

DG121 3Gj2/^^1 D(GL12)/DX1 

DG122 30 ^ 2 / 9 X 2 D(GL12)/DX2 

RAD R Radius 

DRADl 3R/3Xj D(R)/DXl 

0MEGA3 3R/3X2 D(R)/DX2 

Unit 9 Boundary Layer Edge Conditions (7 14.7) 

UlE Ujg XI edge velocity 

U2E U 2 g X2 edge velocity 

PP P Static pressure 

DU1DX2 3Uj/3X2 D(U1E)/DX2 
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DU2DX2 3U2/3X2 D(U2E)/DX2 

DuiDXi au^/aXj d(uie)/dxi 

DU2DX1 aJ2/aXj D(U2E)/DX1 

HR I Rothalpy 

Unit 10 Transformation Variables (7E14.7) 

Y1 or SI XI coordinate or Xl arc length 

Y2 or S2 X2 coordinate or X2 arc length 

( (GAM( I, J) , J=1 ,3) , 1=1 ,3) Direction cosines 
Note: If I0PT9 = 0 Y1 and Y2 are written 
If I0PT9 = 1 SI and S2 are written 
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4.4 Diagnostics 

A number of checks are made during the course of the calculation. If a 
minor error occurs, a DIAGNOSTIC message is printed and the calculation is con- 
tinued. If a serious error occurs, a DIAGNOSTIC message is printed and the cal- 
culation stops. A description of the DIAGNOSTIC message is given below. These 
messages are always of the form 

** DIAGNOSTIC NO. XX ** 

where the number XX refers to one of the errors listed below. 

1) **ERROR** X PT DOES NOT CORRESPOND TO GRID LINE or 
**ERROR** Y PT DOES NOT CORRESPOND TO GRID LINE 

This error is detected in Subroutine INTRID and occurs whenever the computa- 
tional boundary lies outside the boundaries of the input surface coordinates or 
outside the boundaries of the input pressure distribution. This is not a fatal 
error and the calculation continues by extrapolation of the data. However, the 
input data should be examined to be certain it is correct and the computational 
boundaries should be inside the input data boundaries. 

2) **ERROR** 1-D INTERPOLATION ROUTINE CALLED 

HOWEVER NEITHER X NOR Y LIE ON A GRIDLINE 
This error ordinarily should not occur, however, due to round off error, an IF 
test may fail and the error message occurs. This is not a fatal error and the 
calculation continues. 
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3) CONVERGENCE DID NOT OCCUR 

This error is detected in subroutine UITER and occurs whenever convergence of 
the successive substitution iteration did not occur or if only imaginary solu- 
tions exist. This is a fatal error and the calculation stops. Check to see if 
the input pressure distribution is correct and that the input inflow conditions 
are correct. One may also smooth the input data using the I0PT7 option. 
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4.5 Sample Input 

A sample input to the SETUP code is given below card by card as described in 
Section 4.1. To assist the user in setting up a new case comments are inserted 
between blocks of data in lower case letters. These comments of course should not 


be on the input cards . 
Card( s) 0 Title Card 


LSRR76 PRESSURE SURFACE ROTATING RIG (LEADING EDGE START) 

Card(s) 1 Computational Grid Data 
1 40 40 

.48 

Card(s) 2 Computational Boundary Data 
0 21 21 

.01601 0.94699 0.01893 0.92729 

Card( s) 3 Coordinate Surface Data 


3 

50 

7 






0.01000 

0.02000 

0.03000 

0.04000 

0.05000 

0.06000 

0.07000 

0.08000 

0.09000 

0.10000 

0.12500 

0.15000 

0.17500 

0.20000 

0.22500 

0.25000 

0.27500 

0.30000 

0.32500 

0.35000 

0.37500 

0.40000 

0.42500 

0.45000 

0.47500 

0.50000 

0.52500 

0.55000 

0.57500 

0.60000 

0.62500 

0.65000 

0.67500 

0.70000 

0.72500 

0.75000 

0.77500 

0.80000 

0.82500 

0.85000 

0.87500 

0.90000 

0.91000 

0.92000 

0.93000 

0.94000 

0.95000 

0.96000 

0.97000 

0.98000 







0.01892 

0.11828 

0.23656 

0.47311 

0.70967 

0.82795 

0.92730 


-.42037 

-.40957 

-.40300 

-.39908 

-.39723 

-.39723 

-.39909 

-.40302 

-.40933 

-.41594 

-.43076 

-.44336 

-.45392 

-.46258 

-.46945 

-.47461 

-.47813 

-.48006 

-.48045 

-.47932 

-.47670 

-.47261 

-.46705 

-.46003 

-.45157 

-.44165 

-.43026 

-.41739 

-.40301 

-.38708 

-.36958 

-.35047 

-.32973 

-.30730 

-.28313 

-.25717 

-.22931 

-.19947 

-.16755 

-.13341 

-.09689 

-.05780 

-.04139 

-.02451 

-.00716 

.01069 

.02542 

.03137 

.03310 

.03140 







-.43945 

-.42865 

-.42209 

-.41817 

-.41632 

-.41633 

-.41820 

-.42213 

-.42841 

-.43489 

-.44939 

-.46164 

-.47185 

-.48015 

-.48667 

-.49147 

-.49463 

-.49620 

-.49620 

-.49465 

-.49158 

-.48697 

-.48083 

-.47315 

-.46393 

-.45318 

-.44085 

-.42688 

-.41118 

-.39373 

-.37451 

-.35352 

-.33075 

-.30621 

-.27988 

-.25175 

-.22177 

-.18991 

-.15611 

-.12028 

-.08234 

-.04217 

-.02543 

-.00831 

.00920 

.02713 

.04187 

.04782 

.04955 

.04785 







-.46218 

-.45138 

-.44483 

-.44091 

-.43906 

-.43907 

-.44094 

-.44487 
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Card(s) 3 Coordinate Surface Data (Cont'd) 


45108 

-.45739 

-.47142 

-.48321 

-.49295 

-.50078 

-.50682 

-.51116 

51384 

-.51490 

-.51438 

-.51227 

-.50858 

-.50328 

-.49636 

-.48778 

47754 

-.46560 

-.45193 

-.43641 

-.41896 

-.39957 

-.37824 

-.35501 

32990 

-.30297 

-.27425 

-.24375 

-.21150 

-.17748 

-.14170 

-.10411 

06470 

06913 

-.02342 

.06743 

-.00636 

.01101 

.02869 

.04671 

.06145 

.06740 

50768 

-.49688 

-.49031 

-.48639 

-.48453 

-.48453 

-.48639 

-.49032 

49628 

-.50213 

-.51503 

-.52570 

-.53434 

-.54107 

-.54601 

-.54924 

55078 

-.55066 

-.54889 

-.54545 

-.54029 

-.53336 

-.52457 

-.51385 

50108 

-.48617 

-.46906 

-.44971 

-.42815 

-.40445 

-.37870 

-.35101 

32148 

-.29021 

-.25728 

-.22276 

-.18672 

-.14922 

-.11032 

-.07003 

02840 

10830 

.01453 

.10659 

.03207 

.04980 

.06773 

.08587 

.10062 

.10657 

55319 

-.54238 

-.53581 

-.53188 

-.53001 

-.53000 

-.53185 

-.53575 

54132 

-.54659 

-.55810 

-.56744 

-.57477 

-.58021 

-.58387 

-.58578 

58597 

-.58445 

-.58120 

-.57616 

-.56926 

-.56039 

-.54941 

-.53611 

52026 

-.50166 

-.48022 

-.45616 

-.42973 

-.40121 

-.37080 

-.33871 

30503 

-.26993 

-.23348 

-.19578 

-.15691 

-.11692 

-.07589 

-.03386 

00913 

14746 

.05302 

.14576 

.07083 

.08876 

.10683 

.12502 

.13978 

.14573 

57592 

-.56511 

-.55854 

-.55461 

-.55275 

-.55274 

-.55459 

-.55850 

56381 

-.56877 

-.57953 

-.58815 

-.59479 

-.59956 

-.60253 

-.60374 

60322 

-.60096 

-.59691 

-.59103 

-.58320 

-.57331 

-.56115 

-.54648 

52898 

-.50839 

-.48462 

-.45801 

-.42901 

-.39798 

-.36520 

-.33090 

29522 

-.25830 

-.22023 

-.18111 

-.14100 

-.09997 

-.05808 

-.01536 

02813 

16705 

.07237 

.16534 

.09027 

.10827 

.12639 

.14460 

.15936 

.16532 

59501 

-.58420 

-.57764 

-.57371 

-.57185 

-.57185 

-.57371 

-.57763 

58270 

-.58739 

-.59752 

-.60554 

-.61160 

-.61579 

-.61818 

-.61881 

61768 

-.61479 

-.61008 

-.60348 

-.59487 

-.58411 

-.57097 

-.55514 

53625 

-.51396 

-.48822 

-.45947 

-.42829 

-.39515 

-.36036 

-.32422 

28685 

-.24842 

-.20899 

-.16869 

-.12756 

-.08566 

-.04306 

.00021 

04412 

18350 

.08863 

.18179 

.10660 

.12466 

.14282 

.16105 

.17581 

.18176 
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Card(s) 4 Pressure Distribution Data 


2 

10 

7 

4 

1 



0.01600 

0.05600 

0.09100 

0.23800 

0.39300 

0.53600 

0.66300 

0.86200 

0.94700 






0.01892 

0.11828 

0.23656 

0.47311 

0.70967 

0.82795 

0.92730 

0.34170 

0.30740 

0.25390 

0.08700 

0.18180 

-0.08640 

-0.80640 

0.10400 

0.15350 

0.33820 

0.48250 

0.57320 

0.55950 

0.55950 

0.10250 

0.12500 

0.11780 

0.34470 

0.46380 

0.51730 

0.51730 

0.14290 

0.22310 

0.03990 

0.44920 

0.55130 

0.56350 

0.56350 

0.15000 

0.22800 

0.31720 

0.45170 

0.54970 

0.56430 

0.56430 

0.02200 

0.10320 

0.21260 

0.33820 

0.41190 

0.43380 

0.43380 

- 0.12850 

-0.11320 

-0.11840 

0.03430 

0.10320 

0.07000 

0.07000 

- 0.56290 

-0.53210 

-0.47050 

-0.43250 

-0,39840 

-0.44950 

-0.44950 

- 1.17390 

-1.16500 

-1.13420 

-1.10830 

-1.10750 

-1.16750 

-1.16750 

- 2.78660 

-2.73850 

-2.76600 

-2.62450 

-2.69660 

-2.67800 

-2.67800 

Card(s) 5 

Inflow Conditions Data 




4 

6 

7 





0.0000 

0.2819 

0.4950 

0.6669 

0.8149 

0.9480 


0.1000 

0.1000 

0.1000 

0.1000 

0.1000 

0.1000 


0.0200 

0.1182 

0.2366 

0.4731 

0.7097 

0.8280 

0.9273 

0.1000 

0.3500 

0.4500 

0.4500 

0.4500 

0.3500 

0.1000 

Card 6 

Print Options 





4 

0 

1 

1 

31 



Card(s) 7 

Reference 

Flow Variables 




2074.976 

516.035 

107.314 

44.600 

0.52840 

0.00000 


Card(s) 8 

Rotor Parameters 

( load only 

if I0PT8 

= 1) 


107.338 

-1.00000 

0.00000 

0.00000 




3.8039 

4.7158 

0.01829 

0.92733 

0.956 

0.50000 


Card(s) 9 

Upstream Flow Conditions (load 

only if lOPTll = 1) 


0.0200 

2117.000 

519.000 

75.351 

40.000 



0.1182 

2117.000 

519.000 

75.351 

40.000 



0.2366 

2117.000 

519.000 

75.351 

39.000 



0.4731 

2117.000 

519.000 

75.351 

44.600 



0.7097 

2117.000 

519.000 

75.351 

42.000 



0.8280 

2117.000 

519.000 

75.351 

47.000 



0.9273 

2117.000 

519.000 

75.351 

50.000 




0.77100 
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5.0 OPERATION OF BL3D CODE 


5.1 Input Format 


The card input format is described on the pages which follow. In addition 
to the card input, the code also reads input data files contained on units 10, 
11,12 which may be set up by the user or may be set up by the SETUP preprocessor 
code. In the loading of the cards the variable name, column numbers, and descrip- 
tion are given. Integer variables are all right adjusted and real variables may 


have any number of significant figures. 

Card 1 Title card (load 72 alphanumeric characters) 


Card 2 Option Card (612 ,8X ,3F10 ,0) 


IGBAR 


ITURB 


lENRGY 


LINIT 


KINIT 


NITVC 


1-2 


3-4 


5-6 


7-8 


9-10 


11-12 


Levy-Lees option (recommended 0) 

IGBAR = 0 Laminar Levy-Lees variables 

= 1 Turbulent Levy-Lees variables 

Turbulence option 
ITURB = 0 Laminar flow 

= 1 Turbulent/transitional flow 

Wall heat transfer option 
lENRGY = 0 Specify wall temperature 
= 1 Specify wall heat flux 

XI = 0 inflow boundary 

LINIT = 0 Local similarity start 

= 1 Incompressible flat plate start 
= 2 Incompressible stagnation start 
= 3 Input profiles (not Programmed) 

X2 = 0 inflow boundary 
KINIT = 0 Local similarity start 
= 1 Plain of symmetry start 
= 2 Input profiles (not programmed) 

Max. no. iterations (recommended 20) 
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CONV 

21-30 

Convergence tolerance 

(recommended .0001) 

AK 

31-40 

Ratio of adjacent step 

size (recommended 1 

DNl 

41-50 

Step size at wall 

(recommended 0.1) 

Card 

3 Gas Properties 

(5F10.0) 



R 

1-10 

Gas constant 

( air 

1716.0) 

G 

11-20 

Ratio specific heats 

( air 

1.4) 

PR 

21-30 

Prandtl number 

(air 

.72) 

PRT 

31-40 

Turbulent Prandtl number (air 

.92) 

CSTAR 

41-50 

Sutherland's constant 

( air 

198.6) 


Card 4 Flow Conditions (4F10.0) 


XMA 

1-10 

Reference 

Mach number 


TFS 

11-20 

Reference 

static 

temperature 

(deg. R.) 

PFS 

21-30 

Reference 

static 

pressure 

(psf) 

XLREF 

31-40 

Reference 

length 


(ft) 

Card 

5 Virtual orig 

in/heat flux conditions 

(4F10.0) 



XO 

1-10 

XI virtual origin (recommended > 

.0001) 

WO 

11-20 

X2 virtual origin (recommended 

.0000) 

RHWIN 

21-30 

Wall temperature ratio (Tw/Tref) 


QWSIN 

31-40 

Wall heat flux ( ftlb/ft**2/sec) 


Card 6 

Transition points 

(3F10.0) 


XTl 

1-10 

Beginning of transition 


XT2 

11-20 

End of transition 


CISO 

21-30 

Cross flow eddy viscosity parameter 
CISO =1.0 isotropic turbulence 
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Card(s) 7 
LPRO(L) 
KPRO(K) 
KPT 


Print output station numbers (1615) 

Load up to 16 numbers 
Load up to 16 numbers 
Unit 16 print parameter KKPKIWEND 
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The output for 
coordinates, prints 
at selected output 
pages . 

Title Card 

72 alphanumerii 
Echo Input Data 
IE 

IWEND 

IXEND 

KSTOP 

LSTOP 

IGBAR = 0 
= 1 

ITURB = 0 
= 1 

lENRGY = 0 
= 1 

NITCV 


5.2 Output Format 

the BL3D code echoes the input data, lists the computational 
certain calculated terms, and prints the calculated solution 
L,K) points. The output description is given on the following 


Unit 7 Output (file name OUTPUT) 


characters 


Maximum number of points in boundary layer normal direc- 
tion (in X3) 

Maximum number of points in spanwise direction (in X2) 

Maximum number of points in streamwise direction 
(in XI) 

Maximum number of spanwise points at which boundary layer 
calculation is performed 

Maximum number of streamwise points at which boundary 
layer calculation is performed 

Turbulent Levy-Lees variables are used 
Laminar Levy-Lees variables are used 

Turbulent flow (transitional) 

Laminar flow 

Heat flux specified at the wall 
Temperature specified at the wall 

Maximum number of iterations permitted at each station 
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CONV 


Maximum error permitted in F2N at the wall 

XMA 

00 

Freestream Mach number 

TFS 

too 

Freestream static temperature (degrees Rankine) 

PFS 

Poo 

Freestream static pressure (psf) 

R 

R 

Gas constant ( f t**2/sec**2-degree Rankine) 

G 

Y 

Ratio of specific heats 

PR 

Pr 

Prandtl number 

PRT 

^rT 

Turbulent Prandtl number 

CSTAR 

C* 

Constant in Sutherland's viscosity law (degree Rankine) 

AK 


Ratio of adjacent step sizes in boundary layer direction 

DNl 


Step size at the wall in boundary layer 

XLREF 

£ 

Reference length in feet 

XTl 

x^i 

Beginning of transition in XI direction 

XT2 

X^2 

End of transition in XI direction 

CISC 

Co 

Constant used for scaling eddy viscosity in cross flow 
momentum equation for nonisotropic turbulent flows 

DXV 

AXj 

Step size array in streamwise direction 

DWV 

CM 

Step size array in spanwise direction 

XO 

^10 

XI location at first streamwise location 

WO 

^20 

X2 location at first spanwise location 

RHWIN 

Hw/H„ 

Enthalpy ratio at the wall 

QWSIN 

Qw 

Heat flux at the wall 

LPRO 


Controls print in the streamwise direction 
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KPRO Controls print in the spanwise direction 

LINIT Initial conditions; option at L=1 

LINIT = 0 Local similarity 

= 1 Incompressible flat plate 

= 2 Incompressible stagnation line 

= 3 Input profiles 

KINIT Initial conditions option at K=1 

KINIT = 0 Local similarity 

= 1 Plane of S 3 mimetry 

= 2 Input profiles 

Calculated Reference Parameters 

UFS U„ Reference freestream velocity (ft/sec) 

ROFS p„ Reference freestream density (slug/ft**3) 

AMUFS Pgo Reference freestream viscosity (slug/ft/sec) 

REFS Rg„ Reference freestream Reynolds number 

TREF ^ref Reference temperature (deg. Rankine) 

AMURF 1‘ref Reference viscosity ( slug/ f t/sec) 

RE Rg Reference Reynolds number 

CPRIME C' Normalized Sutherland's constant 

Calculated Local Parameters at point (L,K) 

These parameters are referenced to an intrinsic coordinate system composed 
of the boundary layer edge freestream direction and the cross flow (perpendic- 
ular) direction. In addition, certain parameters are referenced to their projec- 
tions onto the Y3 = 0 plane. Parameters in the intrinsic coordinate system are 
labeled (intr.). Parametrs in the Cartesian coordinate system are labeled 
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(Cart.). Length scales are normalized with respect to the reference length XLREF 
and velocities are normalized with respect to the reference freestream velocity 


UFS. 

Xl/L 

Xl 

XI coordinate 

X2/L 

X2 

X2 coordinate 

Yl/L 

Yi 

Y1 Cartesian coordinate or 

Sl/L 

Si 

SI arc length distance 

Y2/L 

^2 

Y2 Cartesian coordinate or 

S2/L 

S2 

S2 arc length distance 

Y3/L 

^3 

Y3 Cartesian coordinate 

CFS 

Cfs 

Streamwise wall friction coefficient (intr.) Zig/Cp^U^^) 

CFC 

Cfc 

Crossflow wall friction coefficient (intr.) 

QW 

Qw 

Wall heat flux 

ST 

S,p 

Stanton number (Q„/(T -T„) 

a W 

TWAD/TFS 


Adiabatic wall temperature ratio 

UTE 

Ule 

Total freestream velocity 

ALPW 

«w 

Limiting streamline flow angle (intr.) 

BETl 

^1 ) 


BET2 

a. 1 

Characteristic directions (intr.) 

CP 

s 

Pressure coefficient 

VIE 

Vie 

VI velocity component (Cart.) 
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V2E 

^2E 

Y2 velocity component (Cart.) 

BETAE 


Boundary layer edge flow direction (Cart 

BETAW 

Sw 

Limiting wall streamline direction (Cart 

YEDGE 

6 

Boundary thickness 

DELS 


Streamwise displacement thickness (intr. 

DELC 


Crossflow displacment thickness (intr.) 

DELUS 

•^us 

Streamwise velocity thickness (intr.) 

DELUC 

*^uc 

Crossflow velocity thickness (intr.) 

THETSS 

®ss 

Streamwise momentum thickness (intr.) 

THETCC 

®cc 

Crossflow momentum thickness (intr.) 

THETSC 

®sc 

AUg momentum thickness (intr.) 

THETCS 

®cs 

Ug AU(, momentum thickness (intr.) 

THETHS 

®HS 

Streamwise enthalpy thickness (intr.) 

THETHC 

®HC 

Crosswise enthalpy thickness (intr.) 

THETDS 

®DS 

Streamwise dissipation thickness (intr.) 

THETDC 

®DC 

Crosswise dissipation thickness (intr.) 

UlE 

"le 

XI component of edge velocity 

U2E 

^2e 

X2 component of edge velocity 

TW/TFS 

T„/T„ 

Wall temperature ratio 

lEDGE 


Index for boundary layer edge 
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Calculated profiles at the point (L,K) 

The following variables are referenced to the intrinsic coordinate system 


and are normalized to 

the local boundary layer edge conditions 

ZETA3 

^3 

Levy-Lees normal coordinate 

X3/L 

X 3 

Normal distance from wall 

US/UTE 

Us/^Te 

Streamwise velocity ratio 

UC/UTE 

Uc/^Te 

Crossflow velocity ratio 

UT/UTE 

Ui/Uie 

Total velocity ratio 

ALPC 

«c 

Flow angle (deg) 

T/TE 

T/Tg 

Static temperature ratio 

VIS 

(p p)/(PgPg) 

Levy-Lees stretching parameter 
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5.3 Input/Output Data Files 

The input data files are units 10, 11, 12 which are read by the main program 
MAIN. The output data files are units 13, 15, and 16. Unit 20 serves as both an 
input and output temporary file. The data on units 10, 11, 12 are writen by the 
SETUP code or the user may setup this data independently. The data on units 13, 
15, 16, and 20 are written by subroutine ROUTPT. A detailed description of these 
files are given below. 

Unit 10 (IDRMI) Metric Data 

The first set of data on this file is the number of points in the computa- 
tional mesh. F0RMAT(3I5) 

IE Number of points in the X3 direction 

IXEND (No. of pts.) - 1 in the XI direction 

IWEND (No. of pts.) - 1 in the X2 direction 

The second set of data is IXEND XI coordinate differences and IWEND X2 co- 
ordinate differences. F0RMAT( 7E14 .7 ) 

DXV(L) ,L=1 , IXEND XI coordinate difference 

DWV(K) ,K=1 , IWEND X2 coordinate difference 

The final set of data consists of the metric data for each mesh point (L,K) 
having a F0RMAT(7E14 .7/7E14 .7) . This data is arrainged to be read by a DO K LOOP 
nested within a DO L LOOP. At each mesh point (L,K) the following data is 
read . 
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H1(K) 

Hi 

XI metric scale coefficient 

H2(K) 

H2 

X2 metric scale coefficient 

DHIDX(K) 

3H^/3Xj 

D(H1)/DX1 

DHIDW(K) 

3H^/3X2 

D(H1)/DX2 

DH2DX(K) 

3H2/3Xj 

D(H2)/DX1 

DH2DW(K) 

3H2/3X2 

D(H2)/DX2 

G12(K) 

Gi2 

Covariant metric tensor component 

GM(K) 

|G| 

Determinant of metric tensor 

DG12DX(K) 

3G^2/3Xj 

D(G12)/DX1 

DG12DW(K) 

3G]^2/^^2 

D(G12)/DX2 

RAD(K) 

R 

Radius 

DRADX(K) 

3R/3Xj 

D(R)/DX1 

DRADW(K) 

3R/3X2 

D(R)/DX2 

0MEGA3(K) 

^3 

Vorticity normal to the surface 

Unit 11 

(IDRM2) Boundary Layer Edge Data 


The data on this file consists of the boundary layer edge conditions for 
each computational point (L,K). This data, like the data above, is read by a DO 
K LOOP nested in a DO L LOOP. At each mesh point (L,K) the following data is 
read. FORMAT(7E14 .7/1E14 .7) . 


U(K) 

Hie 

XI velocity 

component 

V(K) 

H2e 

X2 velocity 

component 

PP(K) 

E 

Static pressure 

DUDW(K) 

3U^g/3X2 

D(U1)/DX2 



V-11 


R85-956834-1 


DVDW(K) 

3U2g/9X2 

D(U2)/DX2 

DUDX(K) 

3Uig/3Xi 

D(U1)/DX1 

DVDX(K) 

9U2g/3X2 

D(U2)/DX2 

HR(K) 


Rothalpy 


Unit 12 (IDRM3) Direction Cosine Data 

The data on this file consists of the direction cosines for each 
computational point (L,K). This data is read by a DO K LOOP nested in a DO L 
LOOP. At each mesh point the following data is read. FORMAT(7E14 .7/5E14 .7) . 


Y1 

Yi 

Y1 

Cartesian coordinate or 

SI 

Si 

Si 

arc length along XI 

Y2 

^2 

Y2 

Cartesian coordinate or 

S2 


S2 

arc length along X2 

Y3 

^3 

Y3 

Cartesian coordinate 


((GAM(I,J,K) ,J=1 ,3) ,1=1 ,3) Direction cosines 
Unit 13 (1DRM4) Crossection Plotting Data 

The following data is output for each mesh point (L,K) by a DO K LOOP nested 
in a DO L LOOP. F0RMAT(2I3 ,9E12 .5) . 


K 


X2 

index 

L 


XI 

index 

Y1(K) 

Yl 

Yl 

Cartesian coordinate or 

Sl(K) 

Si 

XI 

arc length 

VCAP 

Ule 

XI 

edge velocity 

CVCAP 

U2e 

X2 

edge velocity 
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CP 


Static pressure coefficient 

DELS 

«*s 

Streamwise displacement thickness 

DELC 


Crossflow displacement thickness 

THETSS 

®SS 

Streamwise momentum thickness 

THETCC 

®cc 

Crossflow momentum thickness 

ST 

S.p 

Stanton number 

Unit 15 

( IDRM5 ) 

Contour and Vector Plotting Data 


The following data is output for each mesh point (L,K) by a DO K LOOP nested 


in a DO 

L LOOP. 

FORMAT (2 13 

,9E12,5) 

K 


X2 

index 

L 


XI 

index 

Y1 

Yi 

Y1 

Cartesian coordinate or 

SI 

^1 

SI 

arc length along XI 

Y2 

^2 

Y2 

Cartesian coordinate or 

S2 

^2 

S2 

arc length along X2 

Y3 

^3 

Y3 

Cartesian coordinate 

CFS 

Cfs 

Streamwise friction coefficient 

CFG 

Cfc 

Crossflow friction coefficient 

UlE 

Ule 

Y1 

Cartesian component of velocity or 

VIE 

'^le 

XI 

component of velocity 

U2E 

U2e 

Y2 

Cartesian component of velocity or 

V2E 

^2e 

X2 

component of velocity 

EPSIW 

^Iw 

Y1 

Cartesian component of strain or 

F2N(1) 

®lw 

XI 

component of strain 
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EPS2W ^210 Cartesian component of strain or 

G2N(1) ^2 component of strain 

Unit 16 (IDRM6) Special Crossection data 

The following data is output for each mesh point (L,KPT) where KPT is an 

input parameter. F0RMAT(2I3 ,9E12 .5) . 


K 


X2 index 

L 


XI index 

Y1 

SI 

^1 

Y1 Cartesian coordinate or 
SI arc length along XI 

CFS 

Cfs 

Streamwise friction coefficient 

ST 

St 

Stanton number 

TW 

Tw 

Wall temperature 

TWAD 

Ta 

Adiabatic wall temperature 

STE 

Ste 

Local Stanton number 

Unit 20 

(IDRM7) 

Wall Temperature Data 


The folowing data is output/input for each mesh point (L,K) in a DO K LOOP 
nested in a DO L LOOP. If QW = 0, The wall temperature is the adiabatic wall 
temperature and is output onto file 20. If QW > 0 , the adiabatic wall 
temperature is read from file 20. F0RMAT(E12 . 5 ) . 

TW,TWAD Tg Adiabatic wall temperature 
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5.4 Diagnostics 

A number of checks are made during the course of the calculation. If a 
minor error occurs, a DIAGNOSTIC message is printed and the calculation is con- 
tinued. If a serious error occurs, a DIAGNOSTIC message is printed and the 
calculation stops. A description of the DIAGNOSTIC messages is given below. 

1) NEGATIVE TEMPERATURE DETECTED IN MOLVIS AT K= 
where K is the spanwise station where the error is detected. This error is 
detected in subroutine MOLVIS and indicates that a negative temperature is calcu- 
lated in MOLVIS. This negative temperature only occurs if the iteration at the 
(L,K) mesh point does not converge. The usual problem is that negative spanwise 
velocities occur at this spanwise station producing instabilities. The calcula- 
tion is continued by indexing L to the next streanwise station. 

5.5 Sample Input 

A sample input to the BL3D code is given belnw card by card as described in 
Sect. 5.1. To assist the user in setting up a nsw case comments are inserted 
between blocks of data in lower case letters. These comments of course should not 
be on the input cards. 
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Card 1 Title card (load 72 alphanumeric characters) 

GAS TURBINE CASCADE - ROTATING PRESSURE SURFACE 

Card 2 Option Card (6I2,8X,3F10.0) 

1 1 1 0 020 .0001 1.05 0.10 

Card 3 Gas Properties (5F10.0) 

1716.0 1.4 0.72 0.95 198.6 

Card 4 Flow Conditions (4F10.0) 

.0985114 519.0 2117.0 0.9233 

Card 5 Virtual origin/heat flux conditions (4F10.0) 

00.1000 .000 1.0 000.0 

Card 6 Transition points (3F10.0) 

00.1002 00.1003 1.0 

Card(s) 7 Print output station numbers (1615) 

1 5 18 29 39 40 

1 10 20 30 40 
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6.0 DESCRIPTION OF SETUP CODE 
6.1 SETUP Main Program 

The main program for the SETUP code opens and closes all input and output 
files and calls the subroutines which execute the primary tasks depending on the 
input options lOPTl through I0PT12. An overall view of the SETUP code is shown on 
the Global Tree Structure Chart shown on the following pages. This tree struc- 
ture is arranged by tasks and labeled as they appear in the main program SETUP . 
Column 1 shows the main program which calls the subroutines in column 2. Each 
subroutine in column 2 executes a major task and calls the subroutine in column 
3. For any given calculation not all branches of the tree are called. 
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Chart I Global Tree Structure by Task 
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LPL 

KPL 

LBL 

KBL 

LSL 

KSL 

LDL 

KDL 

IDRMl 

IDRM2 

IDRM3 

IDRM4 

IDRM5 

LPLAST 

KPLAST 

LSLAST 

KSLAST 

LBLAST 

KBLAST 

LDLAST 


6.2 SETUP COMMON Block Variables 
PARAMETERS 


(101) 

Max. 

no. Xj mesh points 

(61) 

Max 

no. X 2 mesh points 

(51) 

Max. 

no. input 

Xj boundary points 

(51) 

Max. 

no. input 

X 2 boundary points 

(51) 

Max. 

no. input 

Xj surface points 

(51) 

Max. 

no. input 

X 2 surface points 

(21) 

Max. 

no. input 

X^ inflow points 

(21) 

Max. 

no. input 

X 2 inflow points 


8 

9 

10 
11 
12 


Output unit for metric data 
Output unit for edge conditions 
Output unit for transformation variables 
Output unit for length scales 
Output unit for setup code 
COMMON /INTV/ Integer Variables 
No. input Xj mesh points 
No. input X 2 mesh points 
No. input X^ surface points 
No. input X 2 surface points 
No. input Xj boundary points 
No. input X 2 boundary points 
No. output Xj^ Cp data points 
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KDLAST 

LULAST 

KULAST 

LILAST 

KILAST 

lOPTl 

IOPT2 

IOPT3 

IOPT4 

IOPT5 

IOPT6 

IOPT7 

IOPT8 

I OP T9 

lOPTlO 

lOPTll 

IOPT12 

XI (L) x5; 

X2(K) X2^K 


No. output X 2 Cg data points 
No. input Xj inflow points 
No. input X 2 inflow points 
No. input X Cp data points 
No. input X 2 Cp data points 
COMMON /lOPT/ Input Options 
Surface geometry option 
Grid boundary option 
Pressure data input option 
Mesh distribution option 
Inflow option 
Print option 

Pressure smoothing option 
Rotor/Stator option 
Coordinate output option 
Relative eddy option 
Upstream inflow option 
Not used 


COMMON /INDV/ Independent Variables 
Xj^ coordinate array 
X2 coordinate array 
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DXl(L) 

ax5; 

Xj^step size 

DX2(K) 

^2,K 

X 2 step size 

XIP 

^1 

Xj^ coordinate 

X2P 

X2 

X 2 coordinate 

COMMON /SURFC/ 2D Coordinate Surface Variables 

YY(LK) 

Y^ or Y 2 

Yj^ or Y 2 Cartesian coordinates 

YY3(LK) 

yL 

^3,K 

Y^ Cartesian coordinate 
COMMON /SMIX/ 3D Coordinate Surface Variables 

YIM(L) 


Y]^ Cartesian coordinate of surface 

Y2M(K) 

^2,K 

Y 2 Cartesian coordinate of surface 

Y3M(L,K) 


Y^ Cartesian coordinate of surface 
COMMON /ARCL/ Arc Length Variables 

SSl(K) 

,K 

Xj arc length array 

HSl(K) 

“iK 

Xj^ metric coefficient array 

SLl 

Si 

Xj^ arc length 

SL2 


X 2 arc length 

SSIO 

®10 

Xj^ arc length start 


COMMON /COMB/ Computational Boundary Variables 


Y10X2(K) ,Y20X2(K) Side 1 input boundary 
YlXll(L) ,Y2X11(L) Side 2 input boundary 
Y11X2(K) ,Y21X2(K) Side 3 input boundary 
YIXIO(L) ,Y2X10(L) Side 4 input boundary 
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COMMON /TRANSV/ Transformation Variables 


Y1P,Y2P 

,Y3P Yj,Y2,Y3 

Cartesian coordinates 

DYll 

9Yi/3Xi \ 

First partial derivatives 

mil 

3Y3/9X2 

of Cartesian coordinates 

DY21 

3Y2/3X^ 

with respect to surface 

DY22 

3Y2/3X2 / 

> 

coordinates 

DY31 

3Y3/3X^ I 


DY32 

3 Y 3 / 3 X 2 ^ 


D2Y111 

3^Y3/3Xi3X^ ^ 

i Second partial derivatives 

D2Y112 

92yj/3X33X2 

of Cartesian coordinates 

D2Y122 

32y3/9X23X2 

with respect to surface 

D2Y211 

3^Y2/3Xj3X^ 
2 , 

coordinates 

D2Y212 

3 Y2/3X^3X2 

) 

D2Y222 

3^Y2/3X23X2 I 

i 

D2Y311 

3^Y3/3X39Xj ’ 


D2Y312 

3^Y3/3X33X2 


D2Y322 

3^Y3/3X23X2 



/ 

COMMON /METRIC/ Metric variables 

GL11,GL12,GL13 
GL21,GL22,GL23 
GL31,GL32,GL33 
GUll ,GU12,GU13 

GU21 ,GU22 ,GU23 > G^j Contravariant metric tensor components 
GU31,GU32,GU33 ) 


ij 


Covarient metric tensor components 
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DG 

|G| 

Determinate of metric tensor 

HHl ,HH2,HH3 

H 1 .H 2 .H 3 

Metric scale coefficients 

DHll 

3Hj/9Xi 1 


DH12 

3Hi/3X2 I 

Derivatives of metric scale coeff. 

DH21 

3 H 2 / 3 X 1 1 


DH22 

3H2/9X2 1 


DG121 

3Gi2/9Xi) 

Derivatives of metric tensor 

DG122 

8 Gi 2 / 3^2 1 

COMMON /DIRC/ Direction cosines 

GAM(I,J) 

’'ij 

Direction cosines 

GAMid.J) 


Inverse of direction cosines 


COMMON 

/INPEDG/ Pressure Distributions Variables 

YID(L) 


Output Y coordinate 

Y2D(K) 

'^ 2 ,K 

Output Y coordinate 

Y3D(L,K) 

C \ 
P .X 

Output pressure coeff. 

YII(L) 


Input Y coordinate 

Y2I(K) 

A 

^2 .K 

Input Y coordinate 

Y3I(L,K) 

cL 

S,K 

Input pressure coeff 


COMMON /BLEDGE/ Boundary Edge Variables 

UlE 

Ule 

Xi Edge velocity 

U2E 

U2e 

X 2 Edge velocity 


VI- 7 



R85-956834-1 


pp 

P 

Edge static pressure 

DUIDXI 

8Ule/3Xi ' 


DU1DX2 

3Ui^/3X2 

^ Edge velocity gradients 

DU2DX1 

3U2e/3Xi j 


DU2DX2 

302^/3X2 


CP 

Cp 

Static pressure coefficient 

DCPDXl 

3Cp/3X^ ) 
( 

Static pressure gradients 

DCPDX2 

3Cp/3X2 ) 


ERROR 

e 


DUIDX 

3Uie/3Xi ^ 

^ Explicit velocity gradients 

DU2DX 

3U2e/3Xi ] 

1 

HR 

I 

Edge rothalpy 



COMMON /INFLO/ Inflow Conditions 

U2X2(K) 

^2e,K 

Side 1 input inflow velocity 

U2X1(L) 

U ^ 
^2e 

Side 4 input inflow velocity 

YI2U(K) 

'^2,K 

Side 1 input coordinate 

YIIU(L) 


Side 4 input coordinate 



COMMON /UBLED/ Stored arrays 

U1(K) 

”le,K 

Xj velocity array 

U2(K) 

^2e,K 

X2 velocity array 

CPD(K) 

^P,K 

Pressure coefficient 

ROTHLP(K) 


Rothalphy 
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COMMON /GASP/ Constants 


GASR 

R 

Gas constant 

GAMM 

y 

Ratio of specific heats 

TI 


0.01745329 

PI 


3.1415926 



COMMON /FLOWI/ Input Flow Conditions 

PINFI 

Poo 

Reference static pressure 

TINFI 

t«, 

Reference free stream static temperature 

ROINFI 

Poo 

Reference density 

UINFI 

Uoo 

Reference velocity 

XMINFI 

Moo 

Reference mach number 

BETH 

CO 

Not used 

XLREF 

z 

Reference length 



COMMON /FLOVA^/ Output Flow Conditions 

PINF 

Poo 

Reference static pressure 

TINF 

T 

00 

Reference freestream static temperature 

TREF 

■^ref 

Reference temperature 

HINF 


Reference enthalpy 

PINFL 

Pl 

Rotor inlet static pressure 

HRINFL 

II 

Rotor inlet enthalpy 

DHINF 

AH 

Side 4 delta enthalpy 

ALPII 

“lOO 

Not used 
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PTABSF 

8 

H 

COMMON /FLOWU/ Output Flow Variables 
Absolute total pressure 

HTABSF 

T 

Absolute total enthalpy 

PTABSL 


Absolute rotor inlet total pressure 

VBTI 

VBt 

COMMON /rotor/ Input Rotor Parameters 
Rotor tip velocity 

GOl ,G02,G03 

^oi 

Direction cosines for vorticity vector 

RADHI.RADTI 

> Rip 

Hub/Tip radius 

YRADHI.YRADTI Yjj.Y^ 

Hub, Tip cooordinates 

GAP I 

G 

Rotor midspan gap 

OMEGl ,OMEG2 

, OMEG3 n 

COMMON /ROTORl/ Output Rotor Parameters 
£ Vorticity components in Cartesian coor 

RAD 

R 

Radius 

DRADl 

9R/9Xj ) 


DRAD2 

9R/9X2 1 

Derivat ives 

VBINF 

^B« 

COMMON /ROTOR2/ Rotor Inlet Parameters 
Rotor blade reference velocity 

UP INF 

8 

-e- 

Rotor reference tangential velocity 

VBINFL 

'^BL 

Rotor blade inlet velocity 

UPINFL 

-e- 

Rotor inlet tangential velocity 

SGNR 

(^/-)y 

Sign on radius vector 

SGNU 


Sign on rotor tangential velocity 

DELCP 

ACp 

Rotor CP correction 
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COMMON /UPSTRl/ Input Upstream Flow Variables 

Y2RI(K) 

^2 ,K 

Y 2 input coordinate 

PTRI(K) 

^T,K 

Upstream absolute total pressure 

TTRI(K) 

^T,K 

Upstream absolute total temperature 


COMMON /UPSTR2/ Input Upstream Flow Variables 

UZRI(K) 

^z,K 

Upstream axial velocity 

BTRI(K) 

^1,K 

Upstream relative flow angle 


COMMON 

/UPSTR3/ Output Upstream Flow Variables 

UIA(K) 

^1,K 

Absolute total velocity 

UPA(K) 


Absolute tangential velocity 

UIR(K) 

Ul,K 

Relative total velocity 

UPR(K) 


Relative tangential velocity 

U2S(K) 

^2s,K 

Est. initial rotor spanwise velocity 


COMMON 

/UPSTR4/ Output Upstream Flow Variables 

UZA(K) 

Uz.K 

Absolute axial velocity 

VBB(K) 

'^B.K 

Rotor velocity 

BTA(K) 


Absolute inlet flow angle 

OMS(K) 

^S,K 

Ratio relative vorticity/rotor speed 

RDK(K) 

% 

Rotor radius 


VI-11 



R85-956834-1 



COMMON 

/UPSTR5/ 

Output 

Upstream Flow Variables 

PS(K) 

Pk 

Rotor inlet static pressure 

TS(K) 

Tk 

Rotor inlet static temperature 

HTA(K) 

% 

Absolute 

inlet 

total enthalpy 

PTA(K) 

^T,K 

Absolute 

inlet 

total pressure 

TTA(K) 

^T,K 

Absolute 

inlet 

total temperature 
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BLEDGE 

CUBA 

CUBB 

CUBTWO 

DIRCOS 

EVAL 

INPDAT 

INTRB 

INTRC 

INTRD 

INTREK 

INTREL 

INTRID 

INTR3D 

LAG4PT 

MATRIX 

METCO 

OPEC 

OPTDAT 

PERM 


6.3 Description of SETUP Subroutines 

Integrate boundary layer edge eqs. 

See subroutine CUBTWO 
See subroutine CUBTWO 

Calculate two dimensional cubic spline fit 
Calculate direction cosines and inverse 
Evaluate function and derivatives of spline fit 
Read input data 

Establish the mapping of the boundaries 
Establish surface geometry 
Calculate pressure coeff. and derivatives 
Interpolate inflow velocity on L=1 boundary 
Interpolate inflow velocity on K=1 boundary 
Linear interpolation on Id grid 
Linear interpolation on 2d grid 
Lagrangian 4 point interpolation 
See subroutine CUBTWO 
Calculate metric data 

Calculate one sided proximate end conditions 

Echo input data 

See subroutine CUBTWO 
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REFCON 

RINLET 

ROTOR 

ROTORP 

SETUP 

SMTHD 

SOLVE 

SPLINE 

UITER 


Setup reference conditions 

Calculate rotor inlet conditions 

Estimate upstream total enthalpy and pressure 

Setup pressure distribution for relative eddy approx. 

Main program 

Smooth input CP data 

Compute metric data and b. 1. edge conditions 
Setup derivatives on boundary for surface spline fit 
Iteration algorithm for UlE , U2E 
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Subroutine BLEDGE(L,K) 


Object 

This subroutine integrates the boundary layer edge equations using the input 
pressure (CP) distribution. The boundary layer edge equations are solved at the 
point XI, X2 for the velocities UlE, U2E and the enthalpy HE for the given input 
static pressure. 

Options 


I0PT3 =0 No pressure data 
= 1 CP = CP(X1 ,X2) 

= 2 CP = CP(Y1,Y2) 

= 3 CP = CP(Y1,X2) 

= 4 CP = CP(X1,Y2) 

= 5 Analytical pressure data 

I0PT5 


= 0 

U2E =0.0 

; UlE 

calculated 

from 

= 1 

U2E(0,X2) 

= 0.0 

; U2E(X1,0) 

= 0.0 

= 2 

U2E(0 ,X2) 

= 0.0 

; U2E(X1 ,0) 

input 

= 3 

U2E(0,X2) 

input 

; U2E(X1,0) 

= 0.0 

= 4 

U2E(0,X2) 

input 

; U2E(X1,0) 

input 


List of Symbols 


AO to AlO 
Cl to C6 
D1 to D6 

El to E3 
L, K 


Coefficients of edge equations 
Coefficients of edge equations 
Coefficients of Euler equations 

Coefficients of energy equation 
Index for XI and X2 
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Theory 

The basic equations are described in Sect. 3.11 of Ref 1. This subroutine 
first calculates the coefficients of the differential equations. Then it deter- 
mines if the point L , K is on an inflow boundary (L=l) or (K=l) or an interior 
point . 

If L=1 and I0PT5 = 1 or 2 ; U2E =0 

If L=1 and I0PT5 = 1 or 2 ; U2E is interpolated from the input data using 

subroutine INTREK 


If K=1 and I0PT5 = 3 or 4 ; U2E =0 

If K=1 and I0PT5 = 3 or 4 ; U2E is interpolated from the input data using 

subroutine INTREL 

IF L>1 and K>1 ; U2E and U2E are calculated by integrating the 

boundary layer edge equations 

After U1E,U2E and HE are calculated, all the remaining boundary edge variables 
required by the BL3D code are then calculated. 

References 

1. Anderson, 0. L. : Assessment of a 3-D Boundary Layer Analysis to Predict 
Heat Transfer and Flow Field in a Turbine Passage. UTRC Report R85-956834. 
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Subroutine CUBTWO(Arg. List) 

Object 

This subroutine performs a two dimensional cubic spline fit to unequally 


spaced data. 


Options 


None 


List of Symbols 


N 

Max. X array size 

M 

Max. Y array size 

X(I) 

Input Yl coordinate 

Y(J) 

Input Y2 coordinate 

U(I,J) 

Input Y3 coordinate 

P(I,J) 

DY3/DY1 

Q(I,J) 

DY3/DY2 

R(I,J) 

D2Y3/DY1/DY2 

A(I,J,K,L) 

Spline coefficients 

CDX 

Work arrays 


DY 

ZX 

ZY 

ORD 

UNB 

F 
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Theory 

This subroutine perfortnes a two dimensional cubic spline fit to the input 
blade surface using the theory of Ref. 1. This subroutine is called by subroutine 
SPLINE, and calls subroutines PERM , MATRIX , CUBA , and CUBB. 

References 

1. Spath, H. ; Spline Algorithms for Curves and Surfaces, Utilitas Mathematica 
Publishing Co. Inc., Winnipeg, 1974. 
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Subroutine DIRCOS 

Object 

This subroutine calculates the direction cosine tensor and the inverse of 
the direction cosine tensor for the transformation from the computational to the 
Cartesian coordinates 
Options 
None 

List of Symbols 

See COMMON blocks 
Theory 

The theory for calculating the direction cosine tensor is given in Sect. 
3.10 by Anderson (Ref.l). 

References 

1. Anderson, 0. L. ; Assessment of a 3-D Boundary Layer Analysis to Predict 
Heat Transfer and Flow Field in a Turbine Passage. UTRC Report R85-956834 . 
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Subroutine EVAL(Arg. List) 

Object 

This subroutine evaluates Y3 ,DY3/DY1 .DY3/DY2 ,D2Y3/DY1/DY1 D2Y3/DY2/DY2 
D2Y3/DY1/DY2 at the point (Y1P,Y2P) using the Spline coefficients calculated in 
subroutine SPLINE, 

Options 


None 



List of 

Symbols 



A(I,J,K,L) 

Spline 

coefficients 

LLl 


No, Y1 

input points 

KKl 


No. Y2 

input points 

LSI 


Corner 

point Y1 index 

KSI 


Corner 

point Y2 index 

DY3Y1 

DY3/DY1 

First ( 

derivat ives 

DY3Y2 

DY3/DY2 



DY3Y11 

D2Y3/DY1/DY1 



DY3Y12 

D2Y3/DY1/DY2 

Second 

derivatives 

DY3Y22 

D2Y3/DY2/DY2 



Theory 





The function and its derivatives are evaluated using the formulas derived in 
Ref. 1. For the function we have: 

4 4 

^3 = E E Aijkl (Yi-Yi)NY2-Y2ok)' (1) 

k=l £=1 
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The partial derivatives are evaluated by analytically differentiating Eq. 

( 1 ). 

References 

1. Spath, H.: Spline Algorithms for Curves and Surfaces. Utilitas 

Mathematica Publishing Co., In.c, Winnipeg, 1974. 
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Subroutine INPDAT 


Object 

This subroutine reads all input data for the SETUP code according to the 
input options lOPTl through I0PT12 
Options 


I0PT4 = 0 
I0PT2 = 0 
= 1 
= 2 
= 3 
= 4 


Uniform (X1,X2) computational grid calculated 
Uniform (Y1,Y2) Cartesian grid calculated 
Skewed (Y1,Y2) grid calculated 
One dimensional mapping using sides 1 and 3 
One dimensional mapping using sides 2 and 4 
Two dimensional mapping 


List of Symbols 

See COMMON blocks 
Theory 

This subroutine reads all input data required by the SETUP code and calcu- 
lates certain grids by default since they do not require all the input data for 
general grids. See Options List. 
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Subroutine INTRB 


Object 

This subroutine establishes the mapping of the physical boundaries to the 
computational boundaries. 

Options 

None 

List of Symbols 

See COMMON block variables 
Theory 

The theory for this calculation is given by Anderson (Ref. 1). The physical 
boundaries are known from the input data as ordered pairs of points for the four 
sides of the physical boundary as shown on Fig. 1. These are given by; 

Y1(0,X2P) ,Y2(0,X2P) KBLAST points for Side 1 

Y1(X1P,1) ,Y2(X1P,1) LBLAST points for Side 2 

Yl(l ,X2P) ,Y2(1 ,X2P) KBLAST points for Side 3 

Y1(X1P,0) ,Y2(X1P,0) LBLAST points for Side 4 

where 


( 1 ) 


0.0 < XIP < 1.0 
0.0 < X2P <1.0 


( 2 ) 
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Then considering Side 1, we note that the boundary is only a function of 
X2P. Thus the Side 1 boundary can be written in parametric form as a function of 
X2P by establishing an array X2B(KB) KB=1 ,KBLAST at equal increments 0.0 < X2B < 
1.0 so that there is a one to one correspondence of points. The same procedure 
follows for the other three sides. With the boundaries established in parametric 
form, the Cartesian coordinates (Y1,Y2) of the of the boundary is known as a 
function of the computational coordinates (X1P,X2P). Given the computational 
coordinate (X1P,X2P), the corresponding Cartesian coordinates of the boundary on 
the four sides is obtained by interpolation using the four point Lagrangian 
interpolation formulas contained in subroutine LAG4PT. 

The interior Cartesian coordinate point (Y1P,Y2P) corresponding to the 
computational coordinate point (X1P,X2P) is obtained by the transfinite mapping 
of Gorden and Thiel (Ref. 2) as shown in Sect. 3.10 of Ref. 1. With the mapping 
established, the derivatives of Y1 and Y2 with respect to XI and X2 are obtained 
by analytically differentiating the algebraic mapping formula. Subroutine LAG4PT 
also returns the first and second derivatives of the Cartesian coordinates Yl and 
Y2 of the boundaries with respect to the appropriate computational coordinate XI 
or X2 so that the first and second partial derivatives can be obtained. 
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References 

1. Anderson, 0. L. : Assessment of a 3-D Boundary Layer Analysis to Predict 
Heat Transfer and Flow Field in a Turbine Passage. UTRC Report R85-956834. 

2. Gorden, J. W. and L. C. Thiel: Transfinate Mapping and Their Application to 
Grid Generation, Numerical Grid Generation, F. Thompson Ed., Elesevier 
Publishing Co. Inc. (1982). 
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Subroutine INTRC 


Object 

This subroutine establishes the complete surface geometry by calculating the 
first and second partial derivatives of the surface coordinate Y3 with respect to 
Y1 and Y2 . 

Options 


lOPTl = 0 
lOPTl = 1 
lOPTl = 2 
lOPTl = 3 


Y3 = 0.0 
Y3 = Y3(Y1) 

Y3 = Y3(Y2) 

Y3 = Y3(Y1, Y2) 


List of Symbols 

See COMMON blocks 
Theory 

The Cartesian coordinate of the surface Y3 and the partial derivatives of Y3 
with respect to Y1 and Y2 are obtained by interpolation. 


lOPTl = 0 
lOPTl = 1 
= 2 

lOPTl = 3 


Y3 and all derivatives are zero. 

Y3 and its derivatives are obtained by interpolation 
using a four point Lagrangian interpolation con- 
tained in subroutine LAG4PT. 

Y3 and its derivatives are obtained using a surface 
spline interpolation contained in subroutine EVAL. 
The spline coefficients are calculated first when 
L=1 and K=1 by calling subroutine SPLINE. 
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Subroutine INTRD 

Object 

Calculate pressure coefficient CP and derivatives. 
Options 


I0PT3 = 0 

No 

CP data 

I0PT3 = 1 

CP 

= CP(X1,X2) 

I0PT3 = 2 

CP 

= CP(Y1,Y2) 

I0PT3 = 3 

CP 

= CP(Y1 ,X2) 

I0PT3 = 4 

CP 

= CP(X1,Y2) 

I0PT3 = 5 

CP 

Analytical function 

List of symbols 




AX Logical variable (if TRUE XlP lies on grid line) 

AY Logical variable (if TRUE X2P lies on grid line) 

ONED Logical variable (if TRUE XlP or X2P lies on grid line) 

Theory 

The pressure coefficient CP and its derivatives with respect XI and X2 are 
calculated in this subroutine. If Xl or X2 lie on a grid line, then CP and its 
derivatives with respect to the input variables (Xl or Yl) and (X2 or Y2) are 
calculated using a one dimensional interpolation contained in subroutine INTRID. 
If Xl and X2 do not lie on a grid line, then CP and its derivatives with respect 
to the input variables (Xl or Yl) and (X2 or Y2) are calculated using a linear 
two dimensional interpolation contained in subroutine INTR3D. Then depending on 
the option I0PT3, the partial derivatives with respect to the coordinates are 
calculated . 
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Subroutine INTREK 

Object 

This subroutine interpolates for the inflow velocities along the XI =0.0 
boundary . 

Options 

None 

List of Symbols 

See COMMON Blocks 
Theory 

The input arrays are Y12U(K) and U2X2(K) . Then given Y2 , along the XI = 0.0 
boundary, the inflow velocity U2E is interpolated from the input arrays using the 
four point Lagrangian interpolation contained in subroutine LAG4PT. 
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Subroutine INTREL 

Object 

This subroutine interpolates for the inflow velocities along the X2 = 0.0 
boundary . 

Options 

None 

List of Symbols 

See COMMON Blocks 
Theory 

The input arrays are YIIU(L) and U2X1(L). Then given Yl, along the X2 = 0.0 
boundary, the inflow velocity U2E is interpolated from the input arrays using the 
four point Lagrangian interpolation contained in subroutine LAG4PT. 
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Subroutine INTRlD 


Object 

Interpolates between data points which lie on a gridline. 
Options 

None 

List of Symbols 


X 

Y 

F 

L,K 

XP,YP 

FP 

DFDX 
DFDY 
LD,KD 
LL,KK 
AX, AY 


Ordinates 
Abscissa 
Function F(X,Y) 

Indices of Interpolation 
Interpolation Point 
F(XP,YP) 
d(F(XP,YP))/dx 
d(F(XP,YP))/dy 

Dimensions; X(LD) , Y(KD) , F(LD,KD) 

Actual sizes: X(LL) , Y(KK) , F(LL,KK) 

Logical; AX = .true, denotes 
interpolation is along an X gridline. 

AY = .true, denotes interpolation along a Y gridline. 


Theory 

The function F is evaluated at the point (XP,YP) using simple divided 
differences. This routine differs from INTR3D in that interpolation is strictly 
along a gridline. 
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Object 

This subroutine interpolates 
using linear interpolation of a first 
function F. 

Options 

None 

List of Symbols 
X 
Y 
F 

L,K 

XP,YP 

FP 

DFDX 

DFDY 

DFDXY 

LL,KL 


a point in between gridlines, 
product approximation of the 


Ordinates 

Abscissa 

Function F(X,Y) 

Indices of Interpolation 

Point of Interpolation 

F(XP,YP) 

d(F(XP,YP))/dx 

d(F(XP,YP))/dy 

d2(F(XP,YP))/dy2 

Dimensions: X(LL) 

Y(KL) 

F(LL,KL) 


Subroutine INTR3D 

a function F to 
order tensor 
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Theory 

The function F is locally approximated by the first order tensor pro- 
duct; 

F = (ax+b)*(cy+d) (1) 

where a,b,c and d are chosen so that the local approximation agrees with the 
actual function values at the four corners of grid box whose lower lefthand cor- 
ner has indices (L,K). The returned function value and derivatives are evaluated 
directly from F ' . 
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Subroutine LAG4PT 

Object 

This subroutine interpolates the function and the first and second 
derivatives using a four point Lagrangian interpolation. 

Options 


None 



List of Symbols 



X(I) 

Input 

abscissas 

Y(I) 

Input 

ordinates 

XARG 

Output 

abscissa 

YARG 

Output 

ordinate 

DYDX 

Output 

first derivative 

D2YDX 

Output 

second derivative 


Theory 

The four point Lagrangian interpolation formula is: 

Y = (X - Xd+DHX - X(H-2))(X - X(H-3))Y(I) 

(X(I) - X(I+1))(X(I) - X(I+2))(X(I) - x(l+3)) 

(X - X(T))(X - X(I+2))(X - X(H-3))Y(I+1) 

(X(I) - X(I+1))(X(I+1) - X(I+2))(X(I+1) - X(l+3)) 

( 1 ) 

+ (X - X(T))(X - X(I+1))(X - X(l+3))Y(H-2) 

(X(I) - X(I+2))(X(I+1) - X(l+2))(X(l+2) - X(l+3)) 

(X - X(I))(X - X(I+1))(X - X(l+2))Y(l+3) 

- (XU) - x(i^3))(x(im - x(i4-t) - ) 

The first and second derivatives are obtained by analytically differentiating 
Eq. (1). 
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Subroutine METCO 

Object 

This subroutine calculates the covariant metric tensor components GLIJ , the 
contravar iant metric tensor components GUIJ , the metric coefficients HHI , and 
their derivatives. 

Options 


None 



List of Symbols 



See COMMON blocks 



Theory 



Given all the partial derivatives 

of the transformation calculated 

in 

routine INTRC, the metric data is then 

given by: 


"ij = 

dYK dY^ 
dX£ dXj 

(1) 

«1 

■ 

(2) 

Hz 

= /Gz2 

(3) 


= 1 .0 

(4) 

G = 

det(G£j) 

(5) 

G^j - (G^gG^^ - Gj.^Gj^g)/G 

(6) 


where the groups (i,r,l) and (j,s,t) are cyclic permutations. The derivatives 
are obtained by differentiating Eq. (1). 
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Subroutine 0PEC(X,F,FPA,FPB,NP1) 

Object 

This subroutine calculates the one sided proximate end conditions for 
the surface spline interpolation. 

Options 

None 

List of Symbols 

X Abscissas 

F Ordinates 

FPA Left hand derivative end condition 

FPB Right hand derivative end condition 

NPl Number of ordered pairs of data (NP1>4) 

Theory 

This subroutine evaluates the derivative boundary conditions for the 
cubic spline which is fitted to the input coordinates of the blade surface 
using the theory of McCartin (Ref 1.). This subroutine is called by sub- 
routine SPLINE. 

References 

1. McCartin, B. J. : Theory, Computation, and Application of Exponential 
Splines, Courant Mathematics and Computing Laboratories, October 1981. 
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Subroutine OPTDAT 

Object 

This subroutine echo prints the input data. 

Options 

None 

List of S3rmbols 

See COMMON blocks 
Theory 

This subroutine echo prints the input data, prints the reference condi- 
tions calculated in subroutine REFCON, prints the smoothed pressure CP data 
calculated in subroutine SMTHD, and finally prints the rotor inlet conditions 
calculated in subroutine RINLET. 
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Subroutine REFCON 

Object 

This subroutine calculates the reference conditions to normalize all input 
variables . 

Options 

I0PT8 = 0 (Stator) No coriolus forces 

I0PT8 = 1 (Rotor) Load pressure data hub to tip 
I0PT8 = -1 (Rotor) Load pressure data tip to hub 
List of Symbols 

See COMMON blocks 
Theory 

The reference conditions are defined in Ref. 1 as used in the BL3D 
boundary layer code. 

References 

1. Anderson, O.L.: Assessment of a 3-D Boundary Layer Analysis to Predict 

Heat Transfer and Flow Field in a Turbine Passage. UTRC Report R85-956834. 
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Subroutine RINLET 

Object 

This subroutine calculates the rotor inlet conditions. 


Options 


I0PT8 = 0 (stator) Non coriolus forces 


I0PT8 = 1 (Rotor) Load pressure data hub to tip 

I0PT8 = -1 (Rotor) Load pressure data tip to hub 

lOPTll = 0 Uniform upstream conditions 
lOPTll = 1 Nonuniform upstream conditions 
lOPTll = 2 Approx, upstream cond. for rotor 
List of Symbols 

DUPDR Derivative of absolute tangential velocity 

GAP Gap normal to flow direction 

HGT Span 

OMEGZ Axial component of absolute vorticity 

Axial component of relative vorticity 
Steamwise component of relative vorticity 
Average relative vorticity 
Steam function 

Theory 

The total pressure PTAI(K), total temperature TTAI(K), axial velocity 

UZAI(K), and rotor inlet relative angle BTRI(K) are known from the input 

data. Both the relative and absolute velocity components can be calculated 

from the velocity triangles. Knowing the absolute velocities, the static 
enthalpy and static pressure can then be calculated. 
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The relative rotor inlet vorticity and the initial spanwise velocities 
as estimated using the theory of Dring and Joslyn, Ref. 1. The absolute 
vorticity is given approximately by 


1 9 (RU.) 

<t> 

R 9R 


( 1 ) 


and the relative vorticity by 


“Z 


( 2 ) 


The component of vorticity in the streamwise direction is then 




(3) 


An estimate of the initial spanwise velocity on the rotor blade is obtained 
by solving Poisson's equation using the area average vorticity. 


V 


2 


“ ^ZA 


(4) 


A Fourier series solution of Eq. (4) is given by (Ref. 2) as 


^ZA X (g-x) 



(5) 


where K is odd. The estimated spanwise velocity at the leading edge is 
given by 


VI-39 



R85-956834-1 



References 

1. Dring, R. P. and H. D. Joslyn: The Relative Eddy in Axial Turbine Rotor 

Passages. ASME Paper 83-GT-22, 1983. 

2. Carslaw, H. S. and J. C. Jaeger: Conduction of Heat in Solids. 2nd 

Edition, Oxford University Press, 1959. 
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Subroutine ROTOR 

Object 

This subroutine calculates the upstream total pressure, total enthalpy 
(rothalpy) and local vorticity component normal to the surface 
Options 

I0PT8 = 0 (Stator) Non coriolus forces 

I0PT8 = 1 (Rotor) Load pressure data hub to tip 

I0PT8 = -1 (Rotor) Load pressure data tip to hub 

lOPTll = 0 Uniform upstream conditions 
lOPTll = 1 Nonuniform upstream conditions 
lOPTll = 2 Approx, upstream cond. for rotor 
List of Symbols 

See COMMON blocks 
Theory 

The upstream conditions PTABSI ,TTABSI .UZABSI ,BETARI , are linearly interpo- 
lated from the input data. The velocities are then calculated from the velocity 
triangles and the relative total enthalpy and rothalpy are given by. 


= % - uj - v|/2 

(1) 

I = Rj. - v|/2 

(2) 


The component of vorticity normal to the surface is then 



( 3 ) 
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Subroutine SMTHD 

Object 

This subroutine smooths the input CP data 
Options 

I0PT7 =0 Do not smooth data 
10PT7 > 0 Smooth data IOPT7 times 

List of Symbols 

See COMMON blocks 
Theory 

The input data is first interpolated from the input LILAST*KILAST grid to a 
LDL*KDL grid. Then the CP data is smoothed by averaging with its neighboring 

points using the formula 

^ fk*l * ^ 
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Subroutine SOLVE 

Object 

This subroutine calculates the geometric coordinate data and the boundary 
layer edge conditions for use by the UTRC boundary layer code. 

Options 

I0PT6 Print every lOPTb'th station 

I0PT9 = 0 Output coordinates Y1 and Y2 
I0PT9 = 1 Output arc lengths SI and S2 
List of Symbols 

See COMMON blocks 
Theory 

This subroutine first calculates the computational step size in both the XI 
and X2 directions. A double DO LOOP sweeps through the mesh and calls the sub- 
routines as follows: 


INTRB 

Establish 

mapping of boundaries 

INTRC 

Establish 

surface geometry 

METCO 

Calculate 

metric data 

DIRCOS 

Calculate 

direction cosines 

ROTOR 

Calculate 

component 

upstream total enthalpy, total pressure and vorticity 
normal to surface 

INTRD 

Interpolate pressure distribution 

BLEDGE 

Calculate 

b. 1. edge conditions 
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At each L,K mesh point the following data is written on the output files 


IDRMl 

8 

metric data 

IDRM2 

9 

edge conditions 

IDRM3 

10 

transformation variables 

IDRM4 

11 

length scales 

IDRM5 

12 

setup code output 


R85-956834-1 


Subroutine UITER(Arg. List) 


Object 

This subroutine solves finite difference equations for the boundary layer 
edge velocities UlE , U2E using an iteration algorithm based on successive sub- 
stitution. 

Options 

None 

List of Symbols 


AO to A14 
Cl to C6 
B1 to B12 
L , K 


Coefficients of edge equations 
Coefficients of edge equations 
Coefficients of algebraic equations 
Index for XI and X2 


Theory 

This subroutine solves a set of quadratic algebraic equations by successive 
substitution as described in Sec. 3.8 in Ref. 1. 

References 

1. Anderson, 0. L. : Assessment of a 3-D Boundary Layer Analysis to Predict 
Heat Transfer and Flow Field in a Turbine Passage. UTRC Report R85-956834. 


C 
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7.0 DESCRIPTION OF BL3D CODE 
7.1 BL3D Main Program 

The main program for the BL3D code opens and closes all input and output 
files and calls the subroutines which execute the primary tasks depending on the 
input options described in Section 5.1. An overall view of the BL3D code is 
shown on the Global Tree Structure Chart shown on the following pages. This tree 
structure is arranged by tasks and labeled as they appear in the main program 
BL3D. Column 1 shows the main program which calls the subroutines in Column 2. 
Each subroutine in Column 2 executes a major tasks and calls the subroutine in 
Column 3. For any given calculation not all branches of the tree are called. 
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Chart II Global Tree Structure by Task 



1.0 Open Input/Output File 



2.0 Read Input Data 


RINPUT 



3.0 Initialize Parameters and Arrays 



Start Outer Do Loop (L=l, IXEND) 

4.01 Read Metric Data, Edge Conditons, and Direction Cosines 


Stat Inner Do Loop (K=l , IWEND) 

5.01 Set Boundary Layer Edge Conditions 




5.02 Calculate Molecular Viscosity 


MOLVIS 


5.03 Calculate Turbulent Viscosity 


TURBL 


5.04 Calculate Coefficients of P.D.E. 




AXCOEF 


Start Iteration Loop to Calculate Turbulent Viscosity 


MOLVIS 

TURBL 


6.01 Calculate Coefficients of Matrix Equations 



PDCOEF 
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Chart II Global Tree Structure by Task (Cont'd) 


6.02 Set Boundary Conditions on Matrix Equations 



BNDCND 





6.03 Invert 1 

latrix Eqi 

lations 



SOLVE 





6.04 Set Solution Arrays on Inflow Plains 


SYMSET 

6.05 Check for Convergence 

6.06 Calculate and Write Output Data 
ROUTPT 


6.07 Reinitialize Solution for Next Point 



NEXTPT 





7.0 End Inner Do Loop 

8.0 End Outer Do Loop 


9.0 




Close Input/Output Files 


POSTP 
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7.2 BL3D COMMON Block Variables 
PARAMETERS 


NPHI 

(61) 

Max. no. X 2 mesh points 

NXDIM 

( 101 ) 

Max. no. X^ mesh points 

NORM 

( 201 ) 

Max. no. X^ mesh points 

IDRMO 

(7) 

Printed output unit 

IDRMl 

( 10 ) 

Coordinate and metric data input unit no. 

IDRM2 

( 11 ) 

Boundary layer edge data input unit no. 

IDRM3 

( 12 ) 

Direction cosines data input unit no. 

IDRM4 

(13) 

Contour and vector plot output unit no. 

IDRM5 

(15) 

Crossection plot output unit no. 


COMMON /AXCFCB/ 

Coefficients of Partial Differential Eqs . 

AX(I) 

A 

Coefficients of partial differential eqs. 

Q 

q 

Levy - Lees length scale parameter 

QO 

'lo 

Levy - lees length scale parameter at X 2 = ' 

GBAR(K) 

^ref 

Reference turbulent/laminar viscosity ratio 


COMMON /CNTRL/ Integer Parameters 


IE 


No. X 3 

mesh 

points 

for energy eq. 

IM 

(IE-1) 

No. Xg 

mesh 

points 

for momentum eq. 

IXEND 


No. X^ 

mesh 

points 


IWEND 


No. X 2 

mesh 

points 
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ITURB 


Turbulence model input option 

IGBAR 


Levy - Lees input option 

lENRGY 


Wall boundary condition on energy eq 

ITC 


Flag for negative temperature 

IFL 


Pointer for direct access files 

KTC 


K counter for BL3D failure 

KTCl 


KTC - 1 

KSP 


K counter 

KSPl 


KSP - 1 

KSTOP 


Max. no. calculated X 2 mesh points 

NIT 


Iteration counter 

NITCV 


Max. no. iterations 

CONV 


Convergence criteria 

LSTOP 


Max. no. calculated mesh points 

LINIT 


Xj inflow option 

KINIT 


X 2 inflow option 


COMMON 

/EDGECB/ Boundary Layer edge Conditions 

U(K) 

Ule 

Xj edge velocity 

PP(K) 

Pe 

Edge static pressure 

V(K) 

U2e 

X 2 edge velocity 
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DVDW(K) 

3U2e/9X2 



DUDW(K) 

9Uie/3X2 

Velocity gradients 


DUDX(K) 

3U^g/9X^ 

f 


DVDX(K) 

9U2 3 / 3 X 2 



Y1(K) 


Yj Cartesian coordinate 


Y2(K) 

^2 

Y 2 Cartesian coordinate 


Y3(K) 

^3 

Yj Cartesian coordinate 


GAM(I,J,K) 


Direction cosines 


HR(K) 

I 

Rothalpy (enthalpy) 



COMMON 

/EWINPT/ Input Wall Heat Transfer 

Variables 

EWIN(K,L) 


Wall heat flux array 


RHWIN 

H„/H. 

Wall enthalpy ratio 


QWSIN 

% 

Wall heat flux 



COMMON /FLORF/ Reference Conditions Parameters 

XMA 


Reference freestream Mach no. 

TFS 


Reference freestream static 

temperature 

PFS 


Reference freestream static 

pressure 
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UFS 

00 

Reference freestream velocity 

ROFS 

Poo 

Reference freestream density 

AMUFS 

Poo 

Reference freestream viscosity 

TREF 

*"ref 

Reference temperature 

XLREF 

1 

Reference length scale 

AMURF ’ 

Poo 

Reference viscosity 

RE 


Reference Reynolds no. 

REFS 

^e“> 

Reference freestream Reynolds 

R 

R 

Gas constant 

PR 

P 

r 

Prandtl number (laminar) 

PRT 

^rT 

Prandtl number (turbulent) 

G 

Y 

Ratio of specific heats 

CSTAR 

c* 

Constant in Sutherlands law 

R1 


Not used 

D 


Constant for isentropic flow 

HO 

Hjoo 

Reference rothalpy 

CPRIME 

C 

Constant 

VK 


Constant 

CP 


Pressure coefficient 

P 

COMMON /METCB/ Coordinate Data 


H1(K) 

Hi 

X metric scale coefficient 

H2(K) 

H2 

X metric scale coefficient 
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DHIDX(K) 

/ 3X^ 




DHIDW(K) 

9H, / 9X, 

y Derivatives 

of 

metric scale coefficients 

DH2DX(K) 

9 H 2 / 9 X 1 1 




DH2DW(K) 

9H-, / 9X, 

2 2 y 




G12(K) 

G 12 

Covariant metric component 

DG12DX(K) 

> 

9Gj 2/9^1 j 

\ Derivatives 

of 

covariant metric component 

DG12DW(K) 

9Gi2/3X2 ' 

I 



GM(K) 

|G| 

Determinant 

of 

metric tensor 


COMMON /PCOEF/ Matrix Elements 
AU1F,BU1F,CU1F,BU1V, Matrix elements 

BUI G , BUI T , BU2F , BU2V , 

AU2G,BU2G,CU2G,BU2T, 

AEF.BEF.CEF.AEG, 

BEG,CEG,AET,BET, 

CET,ACF,ACG,BEV, 

DUl ,DU2,DE,DC 

COMMON /PPROCB/ Print Control Parameters 
LPRO(L) L Print control array 

KPRO(K) K Print control array 
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RAD(K) 

R 

DRADX(K) 

8R/ 3X 

DRADW(K) 

3R/3X, 

OMEGA 

0 

0MEGA3(K) 



COMMON /ROTAT/ Rotor Variables 
Radius 

Derivatives of Radius 

Rotor rotational speed 

Vorticity component normal to surface 


COMMON /STPSIZ/ Step Size Variables 


XIl(L) 


Levy - Lees coordinate in Xj^ direction 

DXIl(L) 

«i 

Levy - Lees step size in X^ direction 

DXV(L) 

AXj 

Coordinate step size in X^ direction 

XO 

^10 

Virtual origin for Xj^ direction 

WV(K) 

^2,k 

Coordinate in X£ direction 

DWV(K) 

^2.k 

Coordinate step size in X 2 direction 

WO 

^20 

Virtual origin for X 2 direction 

XN(N) 

?3 

Levy - Lees coordinate in X^ direction 

DN(N) 


Levy - Lees step size in X^ direction 

DNl 

^^^3^5=0 

Initial step size in X^ direction 

AK 

?3 

Stretching parameter 


COMMON 

/TURBCB/ Turbulent Flow Variables 

El(K) 

1 + Ej^/p 

Ratio turbulent/laminar viscosity Xj^ momentum 

E2(K) 

1 + E 2 /P 

Ratio turbulent/laminar viscosity X 2 momentum 


eq. 


eq. 
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E3(K) 

1 + 

p Ratio turbulent/ laminar viscosity energy eq 

EM(K) 

E/u 

Ratio eddy/laminar viscosity 

XTl 

Xti 

Start of transition 

XT2 

Xt2 

End of transition 

CISO 

^SO 

Ratio E2/E1 

XTRAN 

X^ 

Transition coordinate length 


COMMON /VELL/ Boundary Layer Edge Parameters 

VCAP 

Ule 

X^ edge velocity 

TCAP 

Te 

Edge static temperature 

CVCAP 

^2e 

X 2 edge velocity 

PCAP 

Pe 

Edge static pressure 

DTCDX 

3Tg/3X^ 1 
/ 

Derivatives of edge temperature 

DTCDW 

3T^/3X2 1 


DVCDX 

/ 

( 

Derivatives of X^ edge velocity 

DVCDW 

3U^g/3X2 ) 


DCCDX 

3U2e/3Xi 1 

1 

Derivatives of X^ edge velocity 

DCCDW 

3U2e/3X2 j 

1 

DPCDX 

3Cp/3Xi j 

/ 

Derivatives of static pressure 

DPCDW 

3Cp/3X2 ) 

1 

VR 

Uref 

Reference edge velocity 

DVRDX 

3Uref/9^1 

1 Derivatives of reference edge velocity 

DVRDW 

3U^ef/3^2 

) 
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HE 

DHEDX 

DHEDW 

Hxe 

3H^e/9Xi j 
3H^g/3X2 j 

Edge rothalpy 

Derivatives of edge rothalpy 

ROCAP 

Pe 

Edge density 

ROMU 

L 

Edge Levy - Lees parameter 

AMUE 

Pe^Pref 

Edge viscosity 

DCCDWX 

3 ^U 2 g/ 3 X^ 3X2 

Limiting derivative at symmetry boundary 
COMMON /Vise/ Dependent Variables 

VIS(N) 

L 

Levy - Lees stretching function 

V2(N) 


Normal velocity ratio 

V3(N) 

^>K-1 

Normal velocity ratio 

F3N(N) 

( 3F / 3?2 ) 

Derivative of Xj velocity ratio 

G3N(N) 

(3G/3C3),|j_1 

Derivative of X 2 velocity ratio 

F2(N) 


X| velocity ratio 

G2(N) 


X 2 velocity ratio 

T2(N) 


Static temperature ratio 

F2N(N) 

(3F/353>Vk 

Derivative of Xj velocity ratio 

G2N(N) 

00/353),^ 

Derivative of X 2 velocity ratio 

F1(N) 

pL-1 

X^ velocity ratio 

G1(N) 

^>K 

X 2 velocity ratio 
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FIN(N) 

(3F/9^3)^’1 

Derivative of Xj velocity ratio 

GIN(N) 

0G/9C3)^"^ 

Derivative of X 2 velocity ratio 

Y(N) 

^3 

Distance from wall 


F3(N) 


Xj velocity ratio 


G3(N) 

^»K-1 

X 2 velocity ratio 


TTl(N) 

pL-1 

h,k 

Total enthalpy ratio 


TT2(N) 

hVk 

Total enthalpy ratio 


TT3(N) 


Total enthalpy ratio 


TTIN(N) 

(9H/9C3)^K^ 

Derivative of total enthalpy 

ratio 

TT2N(N) 

(9H/9C3)Vj^ 

Derivative of total enthalpy 

ratio 

TT3N(N) 

(9H/953)Vk_i 

Derivative of total enthalpy 

ratio 

TT4N(N) 

(9H/953)VkJ;^ 

Derivative of total enthalpy 

ratio 

F4(N) 

^>K+1 

X]^ velocity ratio 


G4(N) 

r L-1 
^>K+1 

X 2 velocity ratio 


V4(N) 

yL~l 

Normal velocity ratio 
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AXCOEF 

BNDCND 

EDGE 

MAIN 

MOLVIS 

NEXTPT 

PDCOEF 

POSTP 

PREP 

RINIT 

RINPUT 

ROUTPT 

SOLVE 

SYMSET 

TURBL 


7.3 Description of BL3D Subroutines 

Calculate coefficients of differential eqs . 
Set boundary conditions for matrix eqs. 

Set boundary layer edge conditions 
Main program 

Calculate Levy-Lees function 

Reinitialize solution arrays for next point 

Calculate coefficients of matrix equations 

Close input/output files 

Open input/output files 

Initialize parameters and arrays 

Read input data 

Calculate and write output data 
Invert matrix equations 
Set solution arrays at starting planes 
Calculate turbulent viscosity 
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Subroutine BNDCND(Arg. List) 


Object 

This subroutine sets the boundary conditions on the matrix equations. 
Options 


KINIT = 0 
KINIT = 1 
lERNGY = 0 
lERNGY = 1 


Local similarity start at X2 = 0.0 
Plane of symmetry start at X2 = 0.0 
Specify wall temperature 
Specify wall heat flux 


List of Symbols 


EIE 

EW 

EWP 

FIE 

GIE 

L , K 


Edge value of total enthalpy function 
Wall value of total enthalpy function 
Wall value of enthalpy derivative function 
Edge value of XI velocity function 
Edge value of X2 velocity function 
Index for XI and X2 respectively 


Theory 

The boundary conditions are given Sect. 3.6 of Ref. 1. At the edge of the 
boundary layer the conditions are: 


F 


e 


1 .0 


( 1 ) 


Hg =1.0 (2) 

G, . (3) 
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References 

1. Anderson, 0. L. : Assessment of a 3-D Boundary Layer Analysis to Predict 
Heat Transfer and Flow Field in a Turbine Passage. UTRC Report R85-956834. 
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Subroutine EDGE(Arg. List) 

Object 

This subroutine sets the boundary layer edge conditions. 

Options 

LINIT = 0 Local similarity start at XI = 0.0 
LINIT = 1 Flat plate similarity start at XI = 0.0 
LINIT = 2 Stagnation point similarity start at XI = 0.0 
LINIT = 3 Input profiles (not programmed) 

KINIT = 0 Local similarity start at X2 = 0.0 
KINIT = 1 Plane of symmetry start at X2 = 0.0 
KINIT = 2 Input profiles (not programmed) 

List of Symbols 

See COMMON blocks 
Theory 

Given ^2^, 1^ > and P^, the remaining edge conditions are given as 


follows 

Vg = Rfi (1) 

HTe = ^e + V/2 (2) 

"e = ^le " 4e ^ (3) 

Te = ^ [H^e ~ ^ 1 / 2 ] (4) 

Pe = Pe/Te (5) 
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where the derivatives of the edge conditions are obtained by differentiating 
Eqs. 1 through 5. The reference conditions are set as follows: 


Uref = 


Uref = U 2 e 


KINIT#! 


KINIT=1 , K=1 


( 6 ) 

(7) 


The edge viscosity is obtained from Sutherland's law. 

3/2 1 + C 


- T. 


T^ + C 
e 


( 8 ) 
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Subroutine MOLVIS(L,X) 

Object 

This subroutine calculates the Levy-Lees variable and the static 
temperature. 

Options 

None 

List of Variables 

9 

UTESQ Ug Velocity squared 

X Not used 

Theory 

The total velocity is given by: 

^e = (^le^)^ ^^ref^^^ + ^CUi^F) (Uj.gfG) Gi2/(HiH2) (1) 

Then the static temperature is given by 

H = ^ ^ - u2/2] (2) 

and the Levy-Lees function is calcuated from Sutherland's viscosity law where 

the static pressure is constant in the boundary layer 

T - m_ = ::l/2 1 + g 

Pe^^e “ H + 

A check is made on the static temperature during the iteration. 

If H > 0.0 then ITC = 0 
If H < 0.0 then ITC = 1 
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Subroutine NEXTPT(K,L,IFLAG) 

Object 

This subroutine reinitializes the solution arrays for the next point. 
Options 



IFLAG = 

0 

The solution has converged 


IFLAG = 

1 

The solution has not converged 


KINIT = 

1 

Plane of symmetry start on X2 ■ 

List 

of Symbols 


KSP 


K 


KTC 


K counter for convergence failure 

KTCl 


KTC 

- 1 

KSPl 


KSP 

- 1 

IFL 


Record counter 

Theory 




This subroutine writes the current solution on Unit 8 and reads the 
appropriate record from Unit 8 to set up the solution arrays at the next point. 
The record counter IFL is reset to IFL=1 each time L is indexed. Let us define 
the following vector arrays each representing one record. 
f(L,K-l)i = ( H(N) ,F(N) ,G(N) ,V(N) 

f(L,K>2 = ( H(N) ,F(N) ,G(N) ,V(N) 

f(L-l,K)3 = ( H(N) ,F(N) ,G(N) ,V(N) 

f(L-l,K+l)^ = ( H(N) ,F(N) ,G(N) ,V(N) 

The Unit contains the following records (see Fig. 1): 
for IFL = 1,K f(L,K) j 

for IFL = K+1,KSP f(L-l,K) \ 



( 1 ) 
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where KSP < IWEND. Hence to update the arrays we have: 


fj = £2 = f(lFL) 



1 < K < KSP 


( 3 ) 


K = 1 


(4) 


After the new arrays are set, the normal derivative of the dependent variables 
are calculated numerically. 
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Subroutine PDCOEF(L,K,PR) 


Object 

This subroutine calculates the coefficients of the matrix equations. 


Options 


None 

List of Symbols 

AU1F,BU1F,CU1F, Coeff. of XI momentum eq. 

BU1V,BU1G,BU1T, 

DUl 


BU2F,BU2V,AU2G, Coeff. of X2 momentum eq. 

BU2G,CU2G,BU2T, 

DU2 

AEF.BEF.CEF, Coeff. of energy eq. 

BEV,AEG,BEG, 

CEG.AET.BET 

CET.DE 

ACF,ACG,DC Coeff. of continuity eq. 

Theory 

In this subroutine the equations of motion are quasi-linearized and 


discretized to form a set of matrix equations as outlined in Section 3.9 of Ref. 


1 . 


References 


1. Anderson, 0. L. : Assessmnt of a 3-D Boundary Layer Code to Predict Heat 

Transfer and Flow Field in a Turbine Passage. UTRC Report UTRC85-4. 
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Subroutine POST? 

Object 

This subroutine closes all input/output files. 

Options 

None 

List of Symbols 

See COMMON blocks 
Theory 

See Sect. 5.3 for a description of all input/output files for the three 
dimensional boundary layer code. 


VII-22 


R85-956834-1 


Subroutine PREP 

Object 

This subroutine opens all input/output files. 

Options 

None 

List of Symbols 

See CCM40N blocks 
Theory 

See Sect. 5.3 for a description of all input/output files for the three 
dimensional boundary layer code. 
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'ref - 

(3) 

From Sutherland's law we have 



“ref - 2-27 X lO’S + C*) 

(4) 


w„ = 2.27 X 10"® t„/(t„ + C*) 

(5) 

The 

reference total temperature and total enthalpy are 



ttoo = ^ 

(6) 



(7) 

and 

the reference Reynolds numbes are given by 



Re - P.«/Rref 

(8) 


Re« “ 

(9) 


VII-24 



R85-956834-1 


Subroutine RINPUT 

Object 

This subroutine reads input parameters and flow conditions. 
Options 

None 

List of Symbols 

See COMMON blocks 
Theory 

See Sect. 5.1 for input variables used. 
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Subroutine ROUTPT 

Object 

This subroutine calculates output data in an intrinsic coordinate system 
aligned to the freestream flow direction. 

Options 

lERNGY = 0 Specify wall temperature 


lERNGY = 1 

Specify wall heat flux 

List of 

Symbols 


US(N) 

U 

s 

Velocity component tangent to freestream 

UC(N) 

U 

c 

Velocity component normal to freestream 

UN(N) 

^n 

Not used 

UT(N) 

Ut 

Total (magnitude) velocity 

ALPC(N) 

«c 

Angle between 0 and 

ES(N) 


Rate of strain tangent to freestream 

EC(N) 

®c 

Rate of strain normal to freestream 

EN(N) 

®T 

Total (magnitude) rate of strain 

TS(N) 

^s 

Rate of stress tangent to freestream 

TC(N) 

^c 

Rate of stress normal to freestream 

ALP(N) 

ot 

Angle between and nji 

THE 

0 

Angle between n^ and n 2 

UTIE 

^Tle 

Edge velocity tangent to n^ 

UNIE 

%le 

Edge velocity normal to tij^ and n^ 
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ALPE 

“e 

Angle between and 

UTE 

Uxe 

Xotal edge velocity 

UTl 

Uti 

Velocity tangent to nj^ 

UNI 


Velocity normal to nj and ^3 

E 13 

®13 

XI component of (fi 0 03) rate of strain 

E 23 

®23 

X 2 component of (E 0 03) rate of strain 

ETl 

®T 1 

Rate of strain (E 0 ^3) tangent to nj^ 

ENl 

®nl 

Rate of strain (E 0 ^3) normal to and 

ET 


Magnitude of rate of strain (E 0 ^3) 

ALPEE 

“ee 

Angle between (E 0 ^3) and nj^ 

ALPEC 

®ec 

Angle between (E 0 ^3) and 

T 13 

■^13 

XI component of (f 0 ^3) rate of stress 

T 23 

■^23 

X 2 component of (T 0 03) rate of stress 

TAUTl 

Txi 

Rate of stress (f 0 03) tangent to 

TAUNl 

■^Nl 

Rate of stress (X 0 03 ) normal to n^ and 

TAU 

T 

Magnitude of rate of stress 

ALPT 

“x 

Angle between (T 0 ^3) and 

ALPC 

“c 

Angle between (X 0 S3) and 

DY 3 

AY3 

Difference - normal coordinate 
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DELS 

6*/£ 

DELC 

S*/Si 

DELUS 


DELUC 


THETSS 

CO 

CD 

THETCC 


THETSC 

CD 

CO 

O 

THETCS 


THETHS 

®HS'^ ^ 

THETHC 


THETDS 

0Ds/ 

THETDC 


HSS 

’^SS 

YEDGE 

6 

CFS 

Cfs 

CFC 

Cfc 

QW 

Qw 

ST 

St 

TWAD 

Ta 

TW 

T 

w 

RECFAC 

R 

ec 


Streamwise displacement thickness 

Crossflow displacement thickness 

Streamwise velocity thickness 

Crossflow velocity thickness 

Streamwise momentum thickness 

Crossflow momentum thickness 

AU„ momentum thickness 
c s 

U AU momentum thickness 
s c 

Streamwise enthalpy thickness 

Crosswise enthalpy thickness 

Streamwise dissipation thickness 

Crosswise dissipation thickness 

Shape factor 

Edge of boundary layer 

Streamwise wall friction coefficient 

Crossflow wall friction coefficient 

Wall heat flux 

Stanton number 

Adiabatic wall temperature 

Wall temperature 

Recovery factor 
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V1E,V2E 

> 

> 

BETAE 

Pe 

BETAW 

Pw 

EPSIW 

^Iw 

EPS2W 

^ 2 w 

E131 

^13w 

E231 

^23w 

BETl ,BET2 

6 ^ > 62 

Theory 



Cartesian components of edge velocity 
Flow angle at the edge relative to 
Flow angle at the wall relative to Y 2 
Y1 component of wall rate of strain 
Y2 component of wall rate of strain 
E13 at the wall 
E23 at the wall 
Characteristic directions 


Consider a vector J in a nonorthogonal coordinate system (X1,X2,X3) with the 
corresponding unit vectors (n^^.n^jn^). Then, neglecting the normal component, we 
have for the contravariant components: 


A = A^rij + A^?l 2 + 0 n-; 


( 1 ) 


1 2 

The angle between A and A is given by 


cose = Gj 2 /(HjH 2 ) ( 2 ) 

Let use now resolve the vector in components tangent to normal to Aj . 

A = A,j.n 2 + X n-j) + 0 n (3) 

Clearly as shown n Fig. 1 we have: 


Aj = A^ + A^ COS0 


Ajj = (A^ + sine) 


(4) 

(5) 
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|A| = (a| + (6) 

a = tan (Ajj/A^) (7) 

Let us now define an intrinsic coordinate system referenced to the freestream 


direction. Then applying the above formulas for (A = U) we have: 

Ujle ^le ^2e 

^Nle " ^2e (9) 

Ug = (U^ + 

Qg ~ (11) 


Then for any velocity in the boundary layer we have: 


Uti = Uj + U 2 COS0 (12) 

Ujqi = U 2 sin0 (13) 

= (uj + (14) 

a = tan(Ujjj/U.j.p (15) 


Hence the angle between the local flow direction and the free stream flow 
direction is: 

Og = a - Og (16) 

and the velocity components tangent and normal to the freestream direction 
are : 
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Ug = U.J, cos (17) 

= U.J, sin (18) 

For the rate of strain components and rate of stress components we form the 
vectors 

A = E o n^ = + 023 ^2 ( 19 ) 

X = T o S3 = Tj^3 + T23 ^2 (20) 

where the components are given by 


<^0 = ■'“e "U 

( 21 ) 

Cf - Cq/R^ Wo^Vref ^ref ^ 

( 22 ) 

®13 ^o^le P'^Pe 9^/353 

(23) 

®23 = U^ef P/Pe 

(24) 

Ti 3 = Cj U^g ei/Eref 9F/353 

(25) 

■^23 " ^ref ^l^^ref 30/3^3 

(26) 
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and apply the same formula. With the flow angles determined, the characteristic 
directions are given by: 


= MAX(a^) 
82 = MIN(a^,) 


(27) 

(28) 


The boundary layer integral parameters are defined relative to the intrinsic 
coordinates and are given by: 


dY3 = 1/C^ p/pg d?g 


(29) 


P u, 




e e 


j. D 

■ -/ 1 ^ 1 


(30) 


(31) 


6ug/«. 


= / tl - ] dY. 


(32) 


6 uc/ *- 


-/ [1 - u^] dY3 

e 


(33) 


hsl >- ■ / i' - < > ■'''3 


(34) 
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= -/ 


^ ^dY 


(35) 


^ ^ > “''3 


■^e e 


(36) 


^ ^ ^’'3 


(37) 


P U= ,Uo^ 2 


e e 


(38) 


P U„ 2 

Src/^ = / Tu! (uj) *^^3 


e e 


(39) 


0DS/^ = / tV^ *^^3 
Pe^e ® 


(40) 


‘^'^3 

e 

The wall friction coefficients are given by 

Cfs - 2 'sw/<0eC> 

^£c = 2 t„/(p^D,2, 


(41) 


(42) 

(43) 


VII-33 



R85-956834-1 


and the Stanton number by 


St ■ 


Q„/(T,-T„) 


where 


O _ Y~J- u iifi — fAJi ') 

Y RgPj. ^Pgbe '^ref ^ ^3"^” 


The projection of the velocity vector on the Y3 plane are: 


'^le ■ ^11 ^le * ^21 ^2( 


'^2e ^12 ^le * ^22 ^2e 


The limiting wall streamline is defined by 


lim ^2e _ ^2w 
Y 3 ^ 0 


hence 


I 


^Iw ■ ^11 ®13w ^21 ^23w 

^2w " ^12 ®13w * ^22 ®23w 

6w = tan-1 (e2^/ei„) 


(44) 


(45) 


(46) 

(47) 

(48) 

(49) 

(50) 

(51) 


I 
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Subroutine SOLVE(Arg. list) 


Object 

Solve a set of three second order and one first order partial differential 
equations forming a block tri-diagonal matrix. 

Options 

lERNGY = 0 Specify wall temperature 

lERNGY = 1 Specify wall heat flux 


List of Symbols 

EIE 

EW 

EWP 

FIE 

GIE 

AU1F,BU1F,CU1F, 

BU1V,BU1G,BU1T, 

DUl 

BU2F,BU2V,AU2G, 

BU2G,CU2G,BU2T, 

DU2 

AEF,BEF,CEF, 

BEV,AEG,BEG, 

CEG,AET,BET, 

CET,DE 

ACF,ACG,DC 

PF, QF,RF,SF,TF, 

PG, QG,RG,SG,TG, 
PE,QE,RE,SE,TE 


Edge value of total enthalpy function 
Wall value of total enthalpy function 
Wall value of enthalpy derivative function 
Edge value of XI velocity function 
Edge value of X2 velocity function 
Coeff. of XI momentun eq. 

Coeff. of X2 momentum eq. 

Coeff. of energy eq. 


Coeff. of continuity eq. 

Coeff. of upper triangular matrix 
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Theory 

The block tri-diagonal matrix is reduced to the upper triangular matrix 
using block factorization in which the continuity equation is eliminated from the 
matrix reduction by direct substitution. The center block is inverted using 
Cramer's rule. The reduction process starts at the outer edge (N = IE) where the 
boundary conditions are known explicity and continue to (N=l) where the deriva- 
tive boundary condition is evaluated. In the backward substitution the coeffic- 
ients for the normal velocity are evaluated from the other coefficients. 
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Subroutine TURBL(K,L,X .ICNTRL) 


Object 

This subroutine calculates the turbulent viscosity and the turbulent Levy- 
Lees parameter. 

Options 


ICNTRL = 0 Not used at the present time 


ITURB 


IGBAR 


= 1 

= 0 Laminar flow 
= 1 Turbulent/transitional flow 
= 0 Laminar Levy-Lees variables 
= 1 Turbulent Levy-Lees variables 


List of Symbols 


BET X 3 

EMO €g/ y 

L , K 

X X^ 

XTRAN Xrj. 


X3 Distance normal to wall 
Eddy viscosity outer layer 
Mesh point indices 
XI coordinate distance 
Intermittancy factor 


Theory 

A detailed description of the turbulence model is given in Sect. 3.7 of Ref. 
1. The basic turbulence model is that developed by Cebeci and Smith Ref. 2. Then 
in Levy-Lees variables we have the distance normal to the wall is given by 
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L 5 d?3 


The outer layer eddy viscosity is then given by 


= X 


[S L (yg/y^ef^^ 


X3 (c) 


where the terms in the brackets are coded as local variables. For the inner 
layer, the eddy viscosity is 


= l2 I#- 


where the mixing length is given by 


l2 [ k 2 /R^ (/2 gi/Hi) X2 ^2 

I =2 (H L yg/yref^ 


and the damping factor by 


D = 1-exp 


//'*^2?i N "w^e’^^e'v I dU 
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The Clauser outer layer constant is usually taken to be 

X = 0.0168 


and the Prandtl constant by 


K = 0.4357 


( 6 ) 


(7) 


The eddy viscosity is then given by 


e 

y 



= fo 


_o V if _k s .2. 


( 8 ) 


where is the intermittancy 


factor derived from Dhawan and Narasimha, 


Ref . 


3. 


Xrj, = { 1-exp [ -4.6513 (• 


X-X 


JIL. 


X>p2'”^'P2 


)']} 


(9) 


The turbulent Levy-Lees parameter is then determined by the eddy viscosity in 
the outer layer 


e 


ref 


1 


1 + 



IGBAR = 0 
IGBAR = 1 


( 10 ) 
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LIST OF SYMBOLS 

a^, Coefficients of boundary edge equations 


n’ n 

Coefficients of boundary layer equations 


Specific heats 


Cf 

Wall friction coefficient 


C 

Pressure coefficient 


P 



-► 

®i 

Unit vectors Cartesian coordinates 


F 

Xi velocity ratio 

(ui/uie) 


Covariant metric tensor components 


G 

X 2 velocity ratio 

(uz/uref) 

hi, Hi 

Metric scale coefficients 

(h-/l) 

h , H 

Static enthalpy 

(h/ui) 

1 Hrp 

Total enthalpy 

(hx/ui) 

H 

Total enthalpy ratio 

(hj/h^g) 

i , I 

Rothalpy 

(i/ui) 

k 

Thermal conductivity 


1 

Reference length 


L 

Levy - Lees parameter 

(pu/pgyg) 
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M 

Mach number 


-► 

Unit vectors surface coordinates 


P'P > 

Static pressure 


P , P 

Total pressure 

(Px/Poo'^i^ 

Pr> PrT 

Prandtl number laminar/turbulent 

(Cpp/k) 

q 

Levy - Lees length scale parameter 


r , R 

Radius 

(r/1) 


Reynolds number 

(p„u„l/p^e 


Gas constant 


Si. Si 

Arc length distance 

(Si/1) 

St 

Stanton number 


t , T 

Static temperature 

^ ^/*^ref ^ 

Ui, Ui 

Velocity components surface coor. 

(“i/uj 


Axial velocity 



Radial velocity 



Tangential velocity 


'^ref* ^ref 

X 2 reference velocity u^g or U 2 g 

^'^ref/'^co) 

Vi, Vi 

Velocity components Cartesian coor. 

(vi/u„) 

V 

Xg reduced velocity 


Vb» 

Rotor velocity 

(vb/ub> 
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yi> 


Pi 


Y 


Y 


-1 

ij 


6 


* 


e 

e 

^ref 

0 


y 

^ref 

P, Q 

"ij 
0 ) , n 


Surface coordinates 

Cartesian coordinates 

Relative inlet flow angle 

Ratio of specific heats 

Direction cosine tensor components 
Inverse tensor components 

Displacement thickness 
Eddy viscosity 

Viscosity ratio turbulent/laminar 
Reference viscosity ratio 
Momentum thickness 
Static temperature ratio 
Molecular viscosity 
Reference molecular viscosity 
Levy - Lees transformed coordinates 
Density 

Stress components 
Rotor speed 


(x./l) 

(Yi/l) 


^*^/*-ref ^ 

u(t^gf) 

(p/poo> 

(o)l/u„) 


VIIl-3 



R85-956834-1 


Subscripts 

e 

h 

L 

t 

T 

00 

Superscripts 


Tensor Notation 

i.j.k 

i5j,k 


Edge conditions 
Hub 

rotor inlet conditions 
Tip 

Total or stagnation conditions 
Freestream or reference conditions 

Absolute rotor reference frame 

Subscripts covariant tensor components 
Superscripts contravariant tensor components 
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